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Abstract 

We perform a comprehensive study of models of dark matter (DM) in a Universe with a non-thermal 
cosmological history, i.e with a phase of pressure-less matter domination before the onset of big-bang nu- 
cleosynethesis (BBN). Such cosmological histories are generically predicted by UV completions that contain 
gravitationally coupled scalar helds (moduli). We classify the different production mechanisms for DM in 
this framework, generalizing previous works by considering a wide range of DM masses/couplings and allow¬ 
ing for DM to be in equilibrium with a “dark” sector. We identify four distinct parametric regimes for the 
production of relic DM, and derive accurate semi-analytic approximations for the DM relic abundance. Our 
results are particularly relevant for supersymmetric theories, in which the standard non-thermally produced 
DM candidates are disfavored by indirect detection constraints. We also comment on experimental signals 
in this framework, focusing on novel effects involving the power spectrum of DM density perturbations. In 
particular, we identify a class of models where the spectrum of DM density perturbations is sensitive to the 
pressure-less matter dominated era before BBN, giving rise to interesting astrophysical signatures to be looked 
for in the future. A worthwhile future direction would be to study well-motivated theoretical models within 
this framework and carry out detailed studies of the pattern of expected experimental signals. 
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1 Motivation and Introduction 


Apart from its existence, the nature and non-gravitational interactions of dark matter (DM) are still very 
uncertain. The most popular class of dark matter models - Weakly Interacting Massive Particles (WIMPs) 
- rely on two key assumptions to reproduce the observed relic abundance. First, WIMPs are assumed to 
annihilate into Standard Model (SM) particles with an electroweak-scale cross section. Second, the universe 
is usually assumed to be radiation dominated between the end of inflation and matter-radiation equality. 
However, there are no clear indications that either of these assumptions are valid. With regards to the 
former, large regions of WIMP parameter space have been ruled out by various direct and indirect detection 
experiments. With regards to the latter, the earliest evidence for a radiation dominated universe arises during 
BBN, which occurs at temperatures of order an MeV. The energy budget of the Universe has not been probed 
for temperatures above that at the time of BBN. Of course, it is still possible that dark matter is a simple 
WIMP, but because of the above reasons it is well-motivated to go beyond the traditional WIMP paradigm, 
both in terms of DM candidates as well as the production mechanisms for DM. 

A well-motivated alternative to the standard “thermal” cosmological history mentioned above is that of a 
non-thermal cosmological history, in which BBN is preceeded by a phase of pressureless matter domination. 
Such a situation is predicted in many top-down theories for new physics e.g. low-energy limits of supergravity 
and string/M-theory compactihcations. These theories, under some very mild assumptions, contain gravita¬ 
tionally coupled scalars called moduli. When the Hubble parameter drops below moduli masses, moduli begin 
coherent oscillations and behave as pressure-less matter, dominating the energy density of the universe until 
the longest-lived one {(j)) decays to reheat the universe. In these cosmological histories, an electroweak-scale 
Wino provides a natural candidate for supersymmetric (SUSY) DM, provided that the modulus dominated 
phase ends at temperatures below a GeV or so [^[^. However, recent FERMI-LAT and HESS observations of 
Galactic Genter photons have placed severe limits on Wino DM [^|^. If the Wino is stable, satisfying these 
constraints in the cosmological histories mentioned above requires a large hierarchy between the modulus and 
gravitino masses [^. Such a hierarchy is quite unnatural for a broad class of models in which moduli stabi¬ 
lization sets the scale of supersymmetry breaking |5j|^ . This conclusion also holds if the lightest superpartner 
is some more general admixture of MSSM particles [^. A simple way to avoid these constraints is to assume 
that the lightest visible sector superpartner, hereafter referred to as the LOSP, is unstable. 

Motivated by the above statements, this work provides a comprehensive study of relic DM production in 
cosmological histories with a late phase of modulus domination. To perform as general an analysis as possible, 
we go beyond the standard WIMP picture by i). allowing for a wide range of DM masses and annihilation 
cross sections and ii). allowing for the possibility that DM is in kinetic equilibrium with some sector other 
than the visible sector. These two assumptions are well motivated in SUSY theories with an unstable LOSP, 
but can also be true in general. If the LOSP decays, DM is not a visible sector particle; a priori there is 
no reason to expect its DM mass or annihilation cross section to be near the electroweak scale. Moreover, 
if the DM resides in a sector that couples weakly to the visible sector, DM could be in kinetic equilibrium 
with a “dark sector” instead of the thermal bath of visible sector particles. As we will see, our results can be 
straightforwardly reduced to that of the single-sector case despite assumption ii). 

The analysis of DM models in this framework can be effectively separated into three questions. 

• How does one classify production mechanisms for relic DM and accurately compute 

• What is the pattern of experimental and observational signals arising within such a framework? 
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• What kind of DM candidates and interactions naturally arise in well-motivated theories? 


In this paper, we focus primarily on the first question by solving the Boltzmann equations for a two- 
sector system with a late phase of modulus domination. A brief overview of this framework, along with 
the Boltzmann equations describing its cosmological evolution, is presented in Section We then classify 
all potential mechanisms for the production of relic DM, and compute QoMh'^ for a large range of DM 
masses and annihilation cross sections. The entire parameter space of these DM models can be classified into 
four different parametric regimes, each with distinct production mechanisms for relic DM. We derive (semi) 
analytic approximations for in these different regimes, and confirm their validity by comparison with 

the numerical solution. This is the main new result of our work, and is discussed in detail in Section 
Readers may go directly to Section 3.5, which contains a self-contained summary of the above results. In 
Sectionj^ we discuss the implications of our results for UV-motivated SUSY theories where the modulus mass 
is of order the gravitino mass, and show that this framework provides several viable alternatives to MSSM 
dark matter. 


In the remainder of the paper, we briefly discuss the latter two questions listed above. In Section we 
discuss potential experimental signatures of the DM models considered here. A significant portion of the 
parameter space predicts free-streaming lengths characteristic of warm dark matter. Furthermore, we identify 
a class of DM models in which the DM power spectrum is sensitive to the linear growth of subhorizon DM 
density perturbations during the modulus dominated era. This can lead to interesting astrophysical signatures, 
such as an abundance of earth-mass (or smaller) DM microhalos which are far denser than their counterparts 
in standard cosmologies |10| . Finally, Section briefly addresses the third question, and describes work done 
in string theory models that have dark sectors. In a companion paper, we will elaborate further on some 
classes of these models. We present our conclusions in Section The appendices contain technical results 
which will be referred to in the main text. 


2 Overview of Two-Sectors - Models and Cosmology 


The framework considered here consists of two sectors: a visible sector containing SM (and perhaps MSSM) 
particles and a dark sector containing the DM. Both the visible and dark sectors are assumed to have sufficient 
interactions such that thermal equilibrium is separately maintained within the two sectors, whose temperatures 
are T and T' respectively. We assume that there exist very weak portal interactions between the two sectors, 
so that T and T' may not be equal to each other. Finally, we assume that the Universe is dominated by the 
coherent oscillations of a modulus field cj) at some time which is much earlier than when BBN occur^ The 
cosmological framework described above is depicted schematically in Figure The results of our work will 


be straightforward to reduce to the single sector case, see the discussion in Section 3.5 


As denoted in Figure the visible sector contains radiation degrees of freedom R, comprised of relativistic 
particles in equilibrium with the SM bath at temperature T. We also track the abundance of an unstable 
WIMP-like particle X which is in equilibrium with the visible sector. X corresponds to the LOSP in the SUSY 
theories discussed in the introduction. The dark sector is assumed to contain a stable DM candidate X', along 
with dark radiation R'. “Dark radiation” refers to dark sector particles which are in thermal equilibrium and 
are relativistic at a given dark sector temperature T'. Henceforth, visible (dark) sector quantities are denoted 


^In general, there could be many moduli present in the early Universe. In this case, cj) should be thought of as the longest-lived 


modulus. DM produced from shorter-lived moduli will be diluted by entropy production 11 . 
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Figure 1: Schematic representation of the Two-sector Framework under consideration. 

using unprimed (primed) variables. For simplicity and concreteness, we assume that no DM asymmetry 
is present, so DM particles and antiparticles need not be separately tracked in the Boltzmann equations. 
Relaxing this assumption is worth exploring in future studies, see for example [^. Finally we make the 
assumption that Mx, Mx' ^ which is naturally expected for the supersymmetric theories discussed in 
the introduction and in Section |4l 

Before moving on to study the cosmological evolution of this system, it is worth mentioning that there are 
constraints on hidden sector relativistic degrees of freedom during BBN and during recombination, through 
their contribution to the expansion rate of the Universe. These constraints are typically presented in terms 
of the number of effective extra neutrino species AAfgff, which is related to the number of relativistic hidden 
sector degrees of freedom by: 

^Nes{TBBN) = 0.57 g'^{T'QBN)^‘^{TBBN), ^Nes{TcMB) = 2.2 ^‘^{Tcmb) (1) 

where Tbbn ~ 1 MeV, Tqmb ~ 1 eV and ^{T) = (T'/T)^. The current 95% CL bounds are A NesiTBBN) < 
1.44 and AAeff (Lcmb) < 0.4 [^. We will discuss the implications of these constraints for the two sector 
models considered here in Section [3.1.21 

2.1 Cosmological Evolution 

The cosmology of the framework can be studied by writing down the Boltzmann equations for the time 
evolution of the relevant quantities which comprise the total energy density of the Universe. This includes 
the modulus energy density p^j,, the energy density arising from X and X' with number densities nx and 
nx' respectively, and the energy densities of radiation in the visible and dark sector, denoted by pn and pR/ 
respectively. The relevant parameters in the Boltzmann equations turn out to be: 

{Trh,'^x-,Mx-,Mx', {(tv), {a vY , Bx , Bx' , g, g*{T), g'^{T')} . (2) 
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Here Fx is the decay width of the unstable X particle, Mx and Mx' denote the masses of X and X' 
respectively, while {(tv) and {crv)' denote the thermally averaged annihilation cross-section of X and X' 
respectively. g*{T) and g'^{T') are the relativistic degrees of freedom in the visible and dark sectors at a 
given temperature T, T'. The quantities Bx and Bx' denote the branching fractions of the modulus to X 
and X' respectiveljj^ Given the assumption Mx, Mx' 'C m^, rj approximately denotes the fraction of the 
energy density from the modulus going to dark radiation, with the remaining fraction (1 — 77) going to visible 
radiation. Finally, following established convention we define T^h in terms of the decay width of the modulus 
F 0 as follows: 


Trh — -\/F<7,Mpi 


45 


A'K'^g^{TRH) 


1/4 


(3) 


where Mpi = 1.22 x 10^^ GeV is the Planck scale, and gi,{TRH) is the number of relativistic degrees of freedom 


in the visible sector at Trr. We will discuss the physical intepretation of Trr in Section 3.1.2 


A priori, the nine parameters in ([^ can vary over a wide range of values, and could affect the computation 
of the DM relic abundance in a variety of ways. However, we will show that for Fx > O( 1 )F 0 , the DM 
production mechanisms only depend on a subset of the parameters in ([^, in particular: 


{Trh, Btot, g,Mx', {crv)', c/*(r), g'^{T')}, 


(4) 


where Btot = Bx + Bx' if X decays to X' and Btot = Bx' if X does not decay to X'. Note that there is 
no dependence on parameters measuring the attributes of the LOSP X- {Tx,Mx,{crv){\ Furthermore, as 
will be discussed in Section the parameters Trh and are completely determined by the masses and 
couplings of the modulus 4>. Thus these parameters are insensitive to the details of the dark sector. In the 
forthcoming analysis, we find it useful to choose benchmark values for the following parameters: 


Benchmark : Trh = lOMeV, Btot = 0.1, = 50TeV, g = 0.1, gi,{T) = 10.75, g'^,{T') = 10.75 (5) 


The theoretical motivation for these benchmark values will be clear from the discussion in Section m With 
these parameters fixed, the DM abundance will depend only on Mx' and {crv)', and we will see that these can 
take a wide range of values for viable DM production mechanisms. As mentioned above, for most of the paper 
we take Fx > F,^, since this is naturally obtained if Fx is not Planck suppressed. In Appendix [B| however, 
we will briefly discuss the case Fx ^ Fj^,. 


The Boltzmann equations which describe this system are a natural generalization of those which are 
applicable to a single sector framework within a modulus dominated Universe, as studied in 14,^. As 
pointed out in these papers, it is more convenient to define dimensionless variables corresponding to the 
energy and number densities and also to convert derivatives with respect to time to those with respect to the 
(dimensionless) scale factor A = ^, with a/ = Thus, following 


14 


15 


we define: 


psA^ A'^ A^ , A'^ , A^ 

4* = R ^ = nxx^, R = PR'TfC^i ^ 

^RH ^RH ^RH ^RH ^ RH 


H = 


^ R + R' Ex'X' + ExX^ 
^— + - 


A 


Trh 


( 6 ) 


^Note that Bx also includes channels in which <!> decays to X through intermediate states; Bx' is similarly defined. 
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Ex ~ {^“x + and Ex' ~ (-^x/ + are the thermally averaged X, X' energies assuming that 

X and X' are in kinetic equilibrium. The Boltzmann equations in terms of these comoving dimensionless 
variables are: 


^ d A 
~dR 
^dA 


H 


■dX 


H 

H 


dA 

dX' 

~dA 

dFt! 

~dA 




c'P A^!'- (1 - B) (1 - ,)4. + Afp, (rf) 




cV^ThhBx 


L-p AAp,r„„A-V2 (Xp,2 - X^) - aV2x (Fx) 

i^4> ^RH 


(7) 


cI^^ThhBx 




Mp.TrhA-^/^ {avY (x' 


- X'^ ) + 




0^2 {1 - B) rj^ + Mpi 

with Cp = Cl = (|p) and 


eq 


2£x- M' 


Trh 


Y^A^/^X{Tx) 


B = 


AV2 


BxEx + Bx'Ex' 


(r? 


X 


( 8 ) 


Xgg and Xgg are related to the X and X' equilibrium number densities via: 


X 


eq 



3 


gTMx^ 

2x‘^ 


K2 



if Mx > T, 


24 ycgC(3)T3 
Trh) 7r2 


if Mx<.T, 


(9) 


where g counts the degrees of freedom of X and c^ = g (3^/4) for bosonic (fermionic) X. X'^ is given by ([^ 
with primed variables replacing unprimed variables. 


Note that we have assumed in Q that X decays to X'; we neglect X' + ...—)• X inverse decays, as the 
dynamics which fix klx' occur when T' < Mx (see Section 3.2) at which point inverse decays are exponentially 
suppressed. The thermally averaged X decay rate is given by: 


(rx) 


KijMxlT) 

gxK2{MxlTy 


(rx) 


Mx:$>T^ Tx 
9x ' 


( 10 ) 


where Tx is the X decay rate in the X rest frame, and Ki and K 2 are modified Bessel functions of the second 
kind. The quantities (T^) and ^7"^ ^ are respectively the thermally averaged partial widths for X —)■ X'R 
and X — )■ X'R'. In the remainder of this work, we focus on the case where all X decay channels yield X' such 
that ([^ is valid; this corresponds to X and X' both being charged under the DM stabilization symmetry. It 
is also possible for X to instead decay directly to visible radiation, as is the case for R-parity violating SUSY 
models. In this case X and X' are essentially decoupled in the Boltzmann equations, which significantly 
simplifies the analysis. In Sectionwe focus on the more complicated case where X decays to X', and discuss 
how relaxing this assumption affects our results. 

The above differential equations are solved subject to the following initial conditions: 


A = l, 4 > = 4)7 




R = 0, R' = 0, X = 0, X' = 0 


( 11 ) 
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These initial conditions are somewhat unphysical as they imply pn = = 0 at ^4 = 1. However, at early 

times the visible and dark radiation energy densities are subdominant, so this approximation is justihed. 
Hj is the initial value of the Hubble parameter which fixes the initial energy density of the modulus field, 
parameterized by <I> 7 . As we will see, in most cases the DM relic abundance is largely insensitive to the initial 
condition 


3 Solution of the Boltzmann Equations and the Dark Matter Abundance 

Given the system of equations Q, it is possible to numerically solve it for various choices of the parameters 
in Q. However, in order to get a good physical intuition of the qualitatively different mechanisms at play, it 
is advisable to study various approximate (semi) analytic solutions which are applicable in different regions 
of the parameter space. We carry out such an exercise in this section. In Appendix 0 we compare our 
approximations to the full numerical analysis and find very good agreement. 


3.1 Useful Approximations 

We now derive useful approximations which allow us to obtain semi-analytic expressions for in Sec¬ 

tion 3.2 To start with, it is worth noting that remains constant until H ~ T,^ to a very good approximation. 
Thus in the following analysis we set H = throughout the period of modulus domination, considerably 

simplifying the Boltzmann equations. Our strategy will be to use physically well-motivated approximations 
to first solve for R, R' and X, and then use these solutions to study the equation for X'. 


3.1.1 Approximate solutions for <1, R and R' 

~ \/2 

Consider first the Boltzmann equation for <1. With H = , it is straightforward to solve for 




<1/ exp 



( 12 ) 


Thus, as expected, <1 remains approximately constant at and only begins to decay vigorously when the 
dimensionless scale factor satisfies A > A*, with 





(13) 


Now consider the equations for R and R'. As can be seen from Q, in addition to the modulus decay term 
these equations contain the X and X' annihilation terms as well as the X decay term. However, it turns out 
that for Mx,Mx' ^ all these terms are quite sub-dominant compared to the modulus decay term. This 
is because if Mx, Mx' the energy densities of X and X' are subdominant to during the modulus 

dominated era; a more detailed argument for this is presented in Appendix [A) Given this approximation, the 
solutions to Q do not depend on the branching fractions of X. Thus the approximate solutions for R and R' 
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can be found readily by integrating the modulus decay term: 


/ \ 1/2 

J (1 - / (1 - ‘^(^0 dA'-, E!{A) « y^R{A) 

3\2/3 /^y/3 


R{A) : 

-Rfinal = i?(^-^00) « (l-r/)(l-Beff)r 


(14) 


CpJ 




b!. 


'final 


\ — r] 




In the second line of (14), i?finai represents the late time solution for R, i.e when the scale factor A^ A^. 


Note that R ss -Rgnai during the radiation dominated era. We have approximated B as 


-Beflf = 


Bx {Mx^ + 3r,)2)i/2 ^ 2 ^ 


X 1/2 


(15) 


where Tq and approximately correspond to the temperatures at which the integrand (14) peaks. These 
temperatures characterize the transition between modulus and radiation domination, and are defined more 


precisely in Section [3.1.2[ To obtain the result above for i?finah we have expanded the function obtained after 
the integration as a series expansion in ^ with ^ <C 1 and kept the leading term. This can be justified by 
taking 4>/ as given by (0, where Hj is the Hubble parameter when the modulus cj) starts dominating the 
energy density of the Universe. Thus, for Hj = with 7 > l|^one hnds ^ <c 1. 


3.1.2 Temperature-scale factor relation and the “maximum” temperature 

The temperature of a system is measured by the radiation energy density, and the relation between the two 
is given in general by: 


T = 




30 \ _ / 30 \ R^/^ 

Jf)) Uv(r)j (A/Tuh) 


(16) 


In a radiation dominated Universe, it is well known that R^^^ = a) remains constant with time, giving 
T oc a~^. However, the situation is different within a modulus dominated Universe since R^/'^ does not remain 
constant with time. It can be shown that at early times when T 3> Trh, $ ~ <I >7 and the temperatures and 


scale factor are related approximately by 14 

T 


V3"55y* V 9*{T) ) ^ ) ’ 


(17) 


where Tmaxj the maximum temperature attained during modulus domination, is given by: 

1/4 


Trr 


Thus, we see that the temperature has a more complicated dependence on the scale factor compared to that in 


radiation domination. Using the fact that Hj = yT^ with 7 ;i> 1, one finds that Tmax ~ From (f^ 


®We expect 7 1 because the modulus dominates the energy density of the universe when ^ H ^ 
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it is straightforward to relate the visible and dark sector temperatures: 


/ r?g,(r) yg^jT) 

\{l-r])gUT')J T ^ \{l-r,)gi{T')J 


(19) 


Combining (18) and (19) gives for the dark sector. As mentioned in Section 2.1 bounds on N^s at both 
Tbbn ~ 1 MeV and Tcmb ~ 1 eV constrain /Tbbn and Tqj^^/Tcmb, which through (19) can be 
mapped into a constraint on g. Comparing ( [T^ with the bound ([^, we see that the resulting constraint 
on g is insensitive to g'^{T') assuming g'^,{T') ^ 0. Taking g*{TBBN) = 10.75 and g*{TcMB) = 3, the AAgfr 
constraints 0 imply g < 0.20 (BBN) and g < 0.06 (CMB). 


In the presence of dark radiation, Tbh as defined in ([^ no longer corresponds to the visible sector 
temperature when 77 = T^, assuming the modulus has completely decayed ($ = 0). Instead, we define the 
temperatures To, as the visible and dark sector temperatures when 77]^^^ = 


77 


<I >=0 


(87r/3)B2 

Mpi 


{PR + PR')^^'^ 


{8n/3y/^ ( PB 

Mpi \^-v) 


Trh{ 1 - T'j, 


( g*{TD) 

\9'AT'b) 


1/4 

Trh 


( 20 ) 


The bounds from Agfi discussed above imply Tr Trr. Hence, for simplicity we will set g*{TRB) = g*{T£)). 
It is also useful to compute the scale factor Ar, which corresponds to the temperature Tr, T^. We compute 


A/) by substituting T = Td and R = T^R nai in (16): 


Ad = 


r(5/3) (3/2)2/3 (1 _ Heff) 


nl/4 




1/3 




1/3 


( 21 ) 


From the definition of A* in (13), we see that A* 


Ad- 


We caution the reader that the definitions of To, and Ad established above are limited in the following 
sense. The above expressions for Td, T'j^ and Ad were derived from 77 = T,^ assuming that the universe 
has reached radiation domination, i.e. <I> = 0 and R = 7?finah R' = -^finai- However, modulus decay is a 
continuous process which occurs when 77 ~ T^jj, but does not have a well-defined start or end point. Upon 
solving the Boltzmann equations, one finds that when 77 = T^, the modulus has not finished decaying and 
the radiation dominated phase has not yet been reached (7? ^ Tg finai )- In the next subsection we will verify 
this fact graphically, utilizing the full numerical solutions for and R (see Figure below). Despite this 
ambiguity, we find Td, T'j^ and Ad to be useful qualitative proxies for the temperature and scale factor at 
which the universe transitions from the modulus dominated to radiation dominated era. 


3.1.3 Approximate solution for X 

Now consider the Boltzmann equation for X. Motivated by earlier statements, we are interested in the case 
where A is a LOSP with weak scale mass and annihilation cross section; thus Xeq will be exponentially 
suppressed for temperatures of a few GeV. In our analysis, we will mostly consider the situation that the 
LOSP X decays before the modulus (typically much before), i.e. Fx > C>(l)r<^. Such a condition can be 
naturally achieved since the modulus decays by Planck-suppressed operators. In Appendix [B| we will briefly 
consider the case where Fx ^ F,^. 
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In the Boltzmann equation for X', the X decay term grows like thus we are interested in the 

solution for X in the low temperature regimes where Xeq can be neglected (this approximation is justified in 
Appendix]^. With this approximation, the Boltzmann equation for X can be written as: 


Cv ^ f _ 

dlogA yXcrit ^crit {o'v) gxT^^ J y^cvit {o'v) j 


( 22 ) 


where Xcr\t is the critical value required for annihilations to be efficient for a given valne of the Hnbble 
parameter. More precisely, Xcrit is given by: 


A^crit — (^x)crit 


A3 

T3 

^RH 


HA^ 




{av) 


RH M^iTrh {crv) 


(23) 


Now, if the processes for depletion of X (the first and second terms on the right hand side of (22)) and the 
production of X (the third term in the right hand side of (22)) are larger than X itself, then these are each 
faster than the Hubble rate and one rapidly reaches a situation where the two processes cancel each other, 
giving rise to what is known as quasi-static equilibrium (QSE) [^. The QSE solution is found by equating 
the right hand side of (22) to zero: 


Aqse = 


TxA^ 


2 TI^h9x {(tv) 


1 + 


Ag\Bxc]!‘^ $ Trh^ {av) 

cJ^^A3 m^MpiV\ 


1/2 


- 1 


(24) 


Given the criteria described above (|24|), QSE occurs when: 
A3 


Aqse + 


Ex 


{cjv) gx 


RH 


> A 


crit 


& 


A3 

{(Tv) 


A/2 


Bx 


A/2 


Mpi ^ 




> A 


crit ' 


(25) 


Upon inspection, one finds that the QSE condition (25) is equivalent to the familiar condition Txjgx > H. 
Thus, we see that as long as Ex > ^xE^, the QSE condition will be satisfied during the modulus dominated 
era such that X « Aqse for Ex > gxH. 


We can gain further insight into the QSE solution for X by rewriting (24) as: 

1/2 


Aqse = 


Aqse 


ExA3 


‘^'^rh9x {crv) 


1 + 


{crv) 

{crv). 


- 1 


q^^Ex m^Mp'^ 


^ pl/2 ^3 

gxbBx- ^ 


4>: 


{av) 


1/2 


; {av) <C {av)^ 
; {av) > {av)^ 


(26) 


Physically, the QSE solution for X occurs when moduli decay into A, and A decay into A', balance one 
another; this explains the dependence of Xqse on 4>. In the above expression, {av)^ is defined as: 


{crv),, = 


1 


f ^fMpim,pTx 


^dx^xj \^ij 
4.48 X 10^^ GeV^ x 


7^6 

^RH 


K 

9xBxJ ^Ad) 


(27) 


10.75 


/ \/10MeVy/ Ex V / 10.75 \ 

UoTevH Trh j VlO-'^GeVy* \gfiTRH)J 
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Figure 2: Plots of the exact solutions for Y — R and X (normalized to their maximum values) as functions of the scale 
factor A. We have taken Hi = Bx ~ 0 . 1 , {av) = 10 ~^ GeV“^ and Px = 10 “® GeV. All other parameters taken to 

their benchmark values ([^. The dashed vertical line represents the scale factor A = Ad defined in ( |21[ ), which characterizes the 
transition between a modulus dominated and a radiation dominated universe. 


Note that for the benchmark choice of parameters in ([^, and Tx not extremely small, (fju)* is quite large 
(compared to a WIMP cross-section ~ 10“^ — 10“^° GeV~^). We expect the same qualitative conclusion as 
long as the portal coupling is not extremely tiny. Thus for supersymmetric models where X is the LOSP, we 
expect {av) <C (uu)*, and hence 6 ~ 1 in the QSE solution for X in the second line of (26). 

Figure]^ shows a plot of the solutions for the values of i? and X (normalized to their maximum values) 
as functions of the scale factor A for the choice of benchmark parameters as in ([^. As can be seen from (12), 
(14) and (26), respectively, the solutions for ^,R and X do not depend on Mx, Mx' or {av)' to a good 
approximation. Moreover the solution for X depends does not depend on {av) for most models of interest in 
which {av) <C {av)^ as we have just discussed above. 


3.2 Classifying Production Mechanisms for Relic Dark Matter 


We now move on to studying the main quantity of interest - the Boltzmann equation for X', whose solution 
will give us the expression for the relic abundance of dark matter X' in terms of a subset of the 

parameters ([^ appearing in the Boltzmann equations. More precisely, the X' relic abundance is given by: 


^Duh 


2 


fYTj) T, „ , 

pr(t,) " 


Mx: 


A,T, 

R{Tf) TaowTRH 




(28) 


In the above expression, T/ is the temperature at any very late time in which the universe has become 
radiation dominated {Tf <C Td) and the X' comoving abundance has become constant. The parameters 
Fnow ~ 2.35 X 10“^^ GeV and ^Rh? k. 4.17 x 10“^ are the present day temperature and radiation relic density. 
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Taking R{Tf ) « -Rfinai and using (16) to relate Af and Tf, (28) can be written as: 


~ L 


-3/4^^(^/) 


Tn. 


nnh^ L=(l-77)(l-5efr)r(5/3) - 


2/3 


(29) 


In order to derive semi-analytic approximations for X'{Tf) and we will solve the Boltzmann equation 

for X' given the approximations stated in the previous sections. In the following we will show that X'{Tf) oc 
<h/, so ^DMh'^ is insensitive to as mentioned above. 


Using the approximate solutions for <I>, R, R' and X in (12), (14) and (26), respectively, we can reduce the 


system of Boltzmann equations in ([^ to a single ordinary differential equation for the evolution of X': 


dX' M^.Trh {av)'A-^/^ 


dA 


H 


Y"/ 2 

^eq ^ 


+ 


H 


+ U^A'qse' 


1/2 

Cl 


axTl 


(30) 


RH 


where Xqse is defined in (26). Note that if X does not decay to X', the Xq^^ term in (30) is absent. Using 
a similar definition for the critical annihilation for X' as was used for X in (23), one can rewrite (30): 


dX' 

d\ogA 

Kr-AA) = 


■ ^2 ■ 



-b 

i^critj 



HA^ 


X: 


eq 


^crit 


+ 




^crit 

HA^/^ 


1/2 o 
Cp 

d72 


m^Mpi 


$ 


(31) 


{(Jv)'TIjj MpxTRHiav)'' 


where Btot = Bx + Bx' if X decays to vB and Btot = Bx' if X does not decay to X'. Just as for the case of 
X, if the processes of depletion of X' (first term on the right hand side of ([M])) and production of X' (second, 
third and fourth terms on the right hand side of (|31[)) are each greater than X' itself, X' will rapidly reach a 


quasi-static equilibrium (QSE) attractor solution such that terms on the right hand side of (31) cancel among 
themselves: 


-^qse(^) — 


1/2 o 

Cp B^ot 


{avY I 


m^Mpi 


A>]+X, 


I 2 
eq 


1/2 


(32) 


Comparing (31) and (32), and using (26) for the QSE solution for X, we see that the QSE conditions hold 
when: 


^QSE > -^crit • 


(33) 


Note that in contrast to X for {Tx) > E^, X' does not necessarily enter QSE during the modulus dominated 
phase. One reason for this is that in contrast to X, which is assumed to be a WIMP, we are exploring a much 
more general set of possibilities for the mass and interactions of the DM particle X'. 

In order to understand better the broad possibilities that could arise for X', it is important to find the 
conditions necessary for QSE to hold at A Ri Ar. If the QSE conditions hold at A rsi Ar, then the positive 
contribution to X' from modulus decay is annihilated away such that X maintains its QSE value. In this 
case, the final X' abundance is insensitiv^ to modulus decay parameters such as and Btot- Conversely 


Note that as discussed below (271, {av) ^ (fru)* for most models where X is a LOSP, for which 6 « 1 from (|26[). Therefore, 


we have used the expression for Xqse with 6 « 1 in (311. 


®Modulo logarithmic sensitivity, as will be discusssed in Section 3.3.1 
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if QSE does not hold at A ~ Ad, Qx'h'^ will be sensitive to contributions from modulus decay, along with 


other sources for X' production during the modulus dominated era. Comparing (31) and (32), we see that 


requiring XQgg(^£)) > XX^{Ad) places a lower bound on {av)'. Keeping the above statements in mind, it 
is useful to define a critical annihilation cross section such that = ^crit sX A = Ad, to delineate the 

various possibilities: 







1 




D) 


r(5/3)^/® 

3/4 


2.35 


9*{Td)v 
10.75 


3/4 


Trh Mpi ’ 


Mx' ^ T'd 
Mx' <C T'd 


(34) 


( 5* {Trh )) Trh 


(35) 


is the 


where we have approximated H aX A = Ad- In the above expressions, cr = \TiA~^(TTH))' ^ 

degrees of freedom of X', and 6 = 1 (3/4) for bosonic (fermionic) X'. The above expressions were obtained 
by taking X'q —)■ 0 in the Mx' /§> T/, case and X'q^^ = X'g in the Mx' ^ T/, case. 

In the following sections, we will classify production mechanisms for X' according to whether or not 
{(Tv)' |r'=T{, > equivalently whether or not X' annihilations are efficient at T/,. To simplifiy the 

following analysis, we will assume that (av) is temperature independent. The generalization of our results to 
temperature dependent (av)' is presented in Appendix [d{ 


3.3 


Efficient Annihilation at T/,: (av)' > (cv)'^ 


If {av}' > (uv)'^, X' tracks its QSE value until Aqgg drops below at A > Ad- When Aggj. drops 

below annihilations are no longer efficient and the comoving X' abundance becomes constant. The 

dynamics of this process, along with the resulting parametrics for fix', depends on whether or not the freeze- 
out temperature for X', T'pQ, is larger than Th. Here T'pQ is the X' freeze-out temperature, which is computed 
assuming a radiation dominated universe (38). If T'pQ > T'p we can neglect A'q in X'qq^ for T ^ Td', in this 
case Aqgg oc and Aggp, drops below when the modulus decays at T' ~ T'p. If instead T'p > T'pQ, 
X' remains in thermal equilibrium during the onset of radiation domination (Agg^. ~ X'q for T' < r'o)- In 
this case the X' relic abundance is determined by the standard freeze-out mechanism. 


3.3.1 Non-relativistic quasi-static equilibrium (QSE^r) 

Eirst consider the case where T'pQ > T'p such that A'q can be neglected for T' > T'p. Assuming {(tv)' > (cru)/, 
X' tracks Aggp, oc until A> Ad, after which the ratio Agg^^/AE^ begins to drop exponentially due to 
the decay of ‘h according to (12). The final X' value is given by AQgj.(Ac), where Ac is determined by solving 


13 
















X'/X' 



Figure 3: Plot of the exact solution of the Boltzmann equations for X' (normalized to its maximum value) as a function of the 
scale factor A corresponding to the QSEnr mechanism. We have taken Hi = {av)' = 10“® GeV“^ and Mx' ~ 10 GeV, 

with all other parameters set to the benchmark values We have also plotted XA^ (311 and Xggg (32 I. The horizontal dashed 
line corresponding to A = Ac is determined by solving the transcendental equation (361 for Ac- 


the transcendental equation: 


^QSe(^c) — -Xentire) ^ 




log 


Ar 


- 

3^ 


Ar 


3/2 


+ log 


{^V)'c 
K? {(Tv)' 

/q\2/3 

r(5/3) 


- log 


{av)' 


- 1 


(36) 


We have defined Ac = Ac^i^^^ and have used the approximation R{Ac) ~ .Rfinai- Taking k k. 2 gives close 
agreement with the full numerical result. We denote the above mechanism for DM production as QSEdx- 


Upon solving (36) for Ac, it is straightforward to compute Q.DMh'^ using (29) with X'{Tj ) = 

'MxA 


n [QSEn 


^1/2 i3/2^^p(_t,-l/2^3/2 


(Mpim,^ (cru)')^/2 


[nnh^ 


(37) 


r(r(|)(i)V3)VU 

-1 

It*. 

o 

1 
_1 

'MxA 

kc/TJ-'Ts/i J 

Mx'Mpi {avY Trh 

T 

_ ^ now _ 


[QRh^ 


In the above, we have made the approximation e ^ 1* 


2 £p 


1/2 


1, see discussion below (12). Also, in the second 


line, we have used (36) to get rid of the exponential factor in the first line. The factor Ac in the numerator 


depends logarithmically on both (erv)' and {ervYc- 
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3.3.2 Standard freezeout during radiation domination h 


Now consider the case where T'pQ < T'p with {av)' > {crv)'p Then, as discussed above, -^qse(^i 


QSEV-^D^ ^eq(^_D 


which implies that X' is in thermal equilibrium at T' ~ Tp and freezes out at some TpQ < T'p when X'q drops 
below XFp The universe is radiation dominated for T' < T^; thus the X' relic abundance is determined by 
the standard thermal freeze-out mechanism. Furthermore, there are two possible sub cases - i) non-relativistic 
freezeout during radiation domination when Tp > M'^ > TpQ, which we denote as -FO™^, and ii) relativistic 
freezeout when Tp > TpQ > Mx', which we denote as FOl^. The relic abundance in the two cases are 
given by (28) with Tj = Tpo, Tj = TpQ and px'{TpQ) determined by the standard freeze-out calculation. 
Specifically, px'{Tpo) ~ where TpQ is defined by n'^q{TpQ) = H/ (av)'. Assuming non-relativistic 

freeze-out, TpQ is given by solving the transcendental equation: 


Mx 

-L ] 




and the resulting relic abundance is given by: 
n [FO 




AV5 

j^l/4 

1 

1/4 

X p 

'Mx'' 


1 - 

f—l 

1 

CO 

1 _ 

_ 9* {Tfo)9'* ) - 


Mx' Mpi {av)'_ 

_ Tnow _ 


If instead Xp < 3, X' freeze-out occurs relativistically, and: 


(38) 


[QRh^]. (39) 


n h^ [FO 


radi 


30C(3)1 

99*{Tfo) 

3/4 

ci 

'MxP 

TT^ 

.(1 - ^) 9'*{T'fo). 


_9*{Tfo)_ 

T 

_ ^ now _ 


[QRh^ 


(40) 


where = g' {3g'/i) for bosons (fermions). 

Note that although the mechanism for DM production discussed here is standard thermal freezeout, the 
relevant parametric region is very different compared to that of usual thermal WIMP freezeout. In particular, 

here TpQ is smaller than Tp fs Trr ~ ~ ^rh for reasonable choices of parameters. 

This implies that Mx' < x'pTp < 10 x Trr. Furthermore, for the cosmological scenarios described in the 
introduction and in Section one expects Trr to be in the range: few MeV < Trr < 100 MeV. Thus, the 
DM in this case is much lighter than a typical electroweak-scale WIMP, even if the underlying mechanism is 
non-relativistic freezeout during radiation domination (FO'^^). On the other hand, DM undergoing relativistic 
thermal freezeout in the dark sector (TO™*^) is qualitatively similar to the case of neutrino decoupling in the 
visible sector. We reiterate that in all other regions of Mx' and {crv)' parameter space, the standard thermal 
freeze-out calculation will not be valid. 


3.4 


Inefficient Annihilation at Tpi (av)' < (uv)'^ 


We now consider the case where {crv)' < {cv)'^ such that X' is not in QSE for T > Tp. In contrast to the 
previous case, the X' relic abundance will be sensitive to both early-time X' production during the modulus 
dominated era and the modulus branching ratio Btot. Because the annihilation rate r(A') ~ (crv)' is 
much smaller than the Hubble parameter for T' > T'p, the term in (30) can be neglected for T' > T'p. 
The Boltzmann equation for X' becomes linear in this limit, and the contributions to Ddm can be separated 
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into two sources: 


^DM h =VLs 


^2 _i_ o h‘^ 

It -|- i idecay '^ • 


(41) 


The first term, Ojecay^^) is the contribution from modulus and X decays. This term can be computed 
by taking H = and integrating the second term in the RHS of (30) to A = Aj ^ A^. Taking 

p(-2cyV3^/^/') ~ 


ex 


1, equation (|29|) gives: 

.2 ^ ^-3/4 


n 


decay 


/i" 


B, 


tot 


Trh 

_ 


Mx 


- [nnh^ 


(42) 


On the other hand, as the name suggests, flann parameterizes contributions to X' production which arise 


from the annihilation term in (30). This has been discussed in |14| in models with a single sector. There are 
two qualitatively different cases regarding the parameterics of riann 

The first case arises when the DM particle X' attains equilibrium at high temperatures (but {(Jv)' is still 
smaller than {av)'^) and freezes out during modulus domination] hence T^iax > '^fo ^ '^D- Here TpQ is the X' 


freeze-out temperature computed assuming a modulus dominated universe (43). Now, one might naively think 


that both non-relativistic and relativistic thermal freezeout may be possible during modulus domination, just 
as they are during radiation domination (see section 3.3.2). However, as noted in 14 , relativistic freeze-out 


cannot occur during modulus domination if {(tv)' oc (T')” with n < 6. To see this, note that the term 
in (31) corresponding to R'R' —)■ X'X' inverse annihilations scales like X'^^^ {av)' /XX^ oc when 


X' is relativistic. Thus if X' decouples from the thermal bath of dark radiation while relativistic at some 
temperature T^ec’ comoving abundance will continue to grow for T' < due to inverse annihilations, 

provided n < 6. In this work we will only consider n < 6; thus for the models considered here, freeze-out 
during modulus domination occurs only if > T'pQ > Tp and Mx’ > TpQ. 

The second case arises when Tpp > never reaches equilibrium) or when > Tpp > Mx' 

{X' decouples while relativistic). In this case, it turns out that the contribution to DM abundance comes 
predominantly from inverse annihilations via R'R' —)■ X'X', as will be seen shortly. 


3.4.1 Non-relativistic freezeout during modulus domination (FO^°°) 


Let us first consider the case where X' reaches chemical equilibrium and then undergoes freeze-out during 
modulus domination (T^ax > '^fo X T'p). From the arguments above, we note that freezeout can only occur 
when DM is non-relativistic, hence we denote this mechanism as The X' freezeout temperature, 

defined as T'pQ such that n^^,{T'pQ) (av) = H{T'pq), is given by solving the following transcendental equation 


for x'jp = 


/ _ Mxi . 


x'p = In 


( 


V2Vl07r3, 




(43) 


where x'p = Note that the above equation, and hence the parameters T'pQ and x'p, are valid only if 


^FO 

^max ^ '^FO X Tp and Mx' > T'pp, i.e. 1 < x'p < ^ ^. 


Then is given by (29) with X'{Tf) = X' (T'pp): 


[FO 


modi 
nr J 


8ri 


Trh 

3 

/ 4 

Xp 

'Mx 


diiF'pQ) 

Mx_ 


Mx' Mpi (av)' 

T 

_ ^ now 


[^Rh^ 


(44) 
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Figure 4: Plot of the exact solution for X' (normalized to its maximum value) as a function of the scale factor A corresponding 
to the FO^°'^ mechanism. We have taken Hi — lO^^P,^, {av)' = 10“®GeV“^, Mx’ = 10 GeV as in Figure]^ but have instead 
chosen Btot = 0 so that the QSE condition is not satisfied, see (34l. All other parameters set to benchmark values § . For 
comparison we have also plotted the comoving X' equilibrium number density X'^q. 


{(Tv)' > {(tv)'q or Mx' < Mq where: 

{av)'Q{Mx>) = 
Modav)') = 


where x'p is the solution of (43). From (43), it can be seen that the condition Xp > 1 is equivalent to 


2e\/l07r3 

dii^Fo) 

3 


3 

g'g*{TRHy^^ 

_2e\/l07r^ 

g'y^Fo) 


Mx' 


MpiT'^prj 

[MpiT^H {av)' T]] 


(45) 


In addition, x'p must be smaller than Mx'/T'p, which puts an additional constraint on the parameters. Thus 
the parameter space for viable is rather limited, as we will show in Section]^ 

3.4.2 Non-relativistic and relativistic inverse annihilation 

Finally, let us consider the situation when one of the conditions in the previous subsection, i.e. > '^fo ^ 
T'p) or Mx' > T'pq, is not satisfied. In this case, X' is instead populated by R'R' -A X'X' inverse annihilations. 
This occurs if X' never reaches equilibrium for T' < or if X' decouples from the thermal bath while 


relativistic. In either case X''^ <C for T' < Mx', allowing us to neglect the term in (30). Integrating 


the first term on the right hand side of (30) from A = Aq = (8/3)^/® to some scale factor A = Af, one get^ 

(46) 


X'{Af) ^ cj/' Mpi {avY Trh-^ dA 

Jaq 


AV^n' 2 


eq 


H 


°Ao corresponds to the scale factor at which T = Tmax, see (18l. 
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Figure 5: Plot of the exact solution of the Boltzmann equations for X' (normalized to its maximum value) as a function of 
the scale factor A corresponding to the /Anr and lAr mechanisms. We have chosen {av)' = 10“^®GeV“^, Mx' = 10 GeV for 
lAnr, and {(tv)' = 10“^® GeV“^, Mx> = 10~'* GeV for lA^, with Btot = 0 and Hi = lO^^P^; all other parameters taken to their 
benchmark values 


While X' is relativistic, the integrand of (46) grows like in the modulus dominated phase {H ~ 
and falls like ^4“^ in the radiation dominated phase {H ~ ^^R/A). Thus if Mx' > X' production occurs 
predominantly when X’ first becomes non-relativistic, while if Mx' < T^)! X' production occurs predominantly 
at the transition between modulus domination and radiation domination. 

In either case the important dynamics for X' production approximately occurs during modulus domination; 
thus taking H ss we can use (17) to rewrite (46) as: 


X'{T') 


192 


(1257r7)i/2 


9'.\n) 


Trh'^ Mpi {av)' 4>/ 

M3 


dx' x‘ 


'11 / 2 


(47) 


where we have defined x' = Mx'/^'. Here T'^ is defined as the temperature at which the integrand of 
f dx'x'^^n'^3 is peaked, and Tj is a temperature chosen such that X'{T') is essentially constant for T' < T'j. 
For relativistic X', the integrand of (46) peaks at T' ss r^/1.75. Thus we will henceforth take Tj- ss T^/1.75, 


though if Mx' ^ T'j^ the integrand of (46), (47) falls rapidly well before Tj. 


The evaluation of the integral in (47) is different in different regimes. If Mx' > we can evaluate (47) 


assuming X' satisfies Maxwell-Boltzmann statistics. The integral in (47) can then be expressed as: 




I ' '11 ' 

dx X n, 


9'^m3 


1 dx'x'^K2{x' 


eq 


dvri 


"0^ 


(48) 


The function x'‘^K 2 {x')‘^ peaks at « 3.6, corresponding to ps 0.28 Mx'- Thus the maximum X' pro¬ 
duction takes place when X' is non-relativistic, justifying our assumption of Maxwell-Boltzmann statistics. 


Finally, if 


Mv 


> x*, then (48) will be exponentially suppressed, in particular by exp(—2Mx'/^m 


J. We de¬ 
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note the above mechanism of DM production via non-relativistic inverse annihilations as /j4nr- We remind the 


reader that for Mx' > (47) is valid if Tpo > ^max or if T'po > Mx' where TpQ is given by ( [4^ ; otherwise 

is determined by non-relativistic freeze-out during modulus domination as described in Section 3.4.1 


What happens when Mx> < In this case, X' production peaks when X' is relativistic at = T[)/1.75, 


and Fermi-Dirac or Bose-Einstein statistics must be taken into account. The integral in (48) can then be 
expressed as: 

C(3)2cg2Mx'6 ^ 1 756 ^(3)2,^2^^,12 


ixL 


dx X 


r^eq 


TT^ 




(49) 


where again = g' [2>g'/4) for bosons (fermions). We denote the above mechanism of DM production via 
relativistic inverse annihilations as lA^. 


Given (|29|) and (|47|)-(|49|), the relic abundance from inverse annihilations can be readily computed: 


f^ann h \IA^^ 


f^ann h \IA^ 


A&g'‘^Xrf‘ 

125i/27rr5/2 l3/4 
32 0^2^(3)2 (1.75)® 

125i/27ri5/2 7^3/4 
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r73/2 
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fT^\ 
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\' 
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/ 


.Tnow 

/ 

- 

M'x] 



T 

now 


i^Rh^] (50) 
[4lRh^] (51) 


where x = Idx'x'^K 2 {x')‘^ and we have taken g^iTu) = g*{TRH) and g'^{Tp) = gl{Tjr) in (51). Of all the 

^niax 

production mechanisms we have studied, the only scenario where the X' relic abundance depends on is 
/Anr in the case where Mx' > Tmax- 

Note that we have assumed above that {crv)' is independent of temperature. For the QSE^r, 
and the processes which determine the DM relic abundance occur when X' is non-relativistic. Thus for 

these mechanisms, a temperature-independent {crv)' is typically a good assumption for s-wave annihilation 
(p-wave annihilations are considered in Appendix]^. However for lA^, the relevant process responsible for the 
DM abundance (inverse annihilation) takes place when X' is relativisticj^ Since lA^ requires {(tv)' < {crv)'^, 
it is expected that (av)' in this case is schematically given by {crv}' = for some heavy mediator scale A 
and positive integer n. The temperature-independent {crv)' case studied here corresponds to n = 0. Another 
well motivated case is n = 2, corresponding to fermionic X' annihilating via a heavy bosonic mediator. We 
consider this possibility in Appendix [D{ and show that the n = 2 case can be recovered from (51) by making 
the replacement {av)' —> 0.17 x 


3.5 Summary of Results 

In this section, we summarize the results of this section for the benefit of the reader. There are four qual¬ 
itatively distinct parametric regimes for in the framework considered. These different regimes are 

summarized in Table {^'^)' is defined in (34), T'p is defined in (20), and Tmax is dehned in (18). The 


quanitites f'pQ and T'pQ are respectively the X' freezeout temperatures during radiation domination (38) 


and modulus domination (43). Here we briefly review the parametrics for in these different regimes. 


and collect the semi-analytic expressions for derived earlier. In the following expressions we will set 


^The process responsible for DM abundance for does take place when X' is relativistic, but in this case the DM 


abundance is independent of (crv)', see (40l. 
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DM Production Mechanism 

Parametric Region 

I. Efficient Annihilation at T'j^ 

{(Tv)' > {av)'^ 

A. Non-Relativistic QSE 

[QSEc^c] 

Mx' > 

B. FO During Radiation Domination 

Mx' < x'pT'j^ 

II. Inefficient Annihilation at 

{av)' < (av)'^ 

FO During Matter Domination -|- 
A.. 

Production from Modulus Decay 

[For"] 

> T'pq > T'p 
[{av)' > (cjQ'o] 

^ Inverse Annihilations(i?'i?' —)• X'X') -|- 
Production from Modulus Decay 
[lAnr&IAr] 

IIA condition not satisfied 
[{av)' < (cjQo] 


Ta ble 1: Summary of t he different parametric regimes for 
in (341 and (crv)^ in (45l. The temperatures is defined 


ao ___ _ 

(181, To in (201, and T^q above (431. 


m 


g*{TjiH) = 9*{Td) = 10.75, which is the SM value for g^{T) at T ~ 10 MeV. We also assume a fermionic DM 
candidate and set g' = 2. Note that the various mechanisms are valid in different parameteric regions; this is 
reflected in the different fiducial values for Mx' and {ov)' chosen in the expressions below. In Appendix 
we compare our approximate expressions with numerical solutions to the Boltzmann equations Q and hnd 
close agreement. 


I.A: Non-Relativistic QSE {QSE-ay): 

DM annihilations are large enough to drive X' to its quasi-static equilibrium (QSE) value until T' is 
close to T'jy, soon after which QSE is lost and the comoving DM abundance becomes constant. The relic 
abundance in this regime is given by (37): 


D [Q5Enr] 


5.2 X (1 - 77)-3/" 


/ Mx' \ /10MeV\ /10-®GeV-2\ 

VlOGevJ V Trh J V (^' J 


(52) 


Ac is dehned below (36), and lies in the range : 1 < Ac < 5 for < (av)' < 10® QSE^^ is 

the precise generalization of the “non-thermal WIMP miracle” studied in 00 . and also captures the 
sub-dominant logarithmic dependence on {av)' and {av)'^ via Ac which was not considered in \l 


• I.B: Preeze-out during radiation domination Sz FOl^^): 

X' tracks its equilibrium value until after T' ~ T^, and freezes-out after the modulus decays and the 
Universe becomes radiation dominated. Both non-relativistic {EO'j^) and relativistic {FOl^) thermal 
freezeout are possible. FO'^^ is the dark sector analogue of standard WIMP freeze-out during radia- 
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tion domination, while is the dark analogue of neutrino decoupling in the visible sector. This 

mechanism occurs only for Mx' ^ T '^; see Table The relic abundances are given by (39) and (40): 


radi 
nr J 


n [FO 

n [FOl^ 


0.13 X 


100 X 


rf 


1/4 


1/4 


Xp 

ITT 


Mx' 

IKeV 


IQ-^ GeV 
{av)' 


-2 


(53) 

(54) 


x'p is defined in (38) and captures the standard logarithmic sensitivity to {(tv)' for thermal freezeout. 


II.A: Freeze-out during modulus domination and production from modulus decay 

X' reaches its equilibrium value and then freezes out during the modulus dominated phase. After freeze¬ 


out, modulus decay continues to populate X' until T < Tp. As discussed in Section 3.4 non-relativistic 
freeze-out during modulus domination occurs only if T^ax > '^'fo ^ '^'d > TpQ. This implies 


I < x'p < -7j^, and (itu)q < {av)' < {av)'^ where {av)'^ is given in (45). The relic abundance is given by 


DM h'^ = f^decay h'^ + f^ann where Udecay and flann are given respectively by (42) and (44) : 
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x'p in (56) is defined in (43) and is logarithmically sensitive to {av)'. 


II.B: Inverse annihilation and production from modulus decay (/Anr&TAr): 

X' does not undergo freezeout during modulus domination. DM production takes place predominantly 
by inverse annihilations as well as production from modulus decay. Specifically, D dm h'^ = f^decay + 


f^ann where Ddecay IS given by (55), while Dannfi^ gets contributions from inverse annihilations. 
There are two different parametrics for Dann h'^ depending on whether Mx' > T'p or vice versa. 

(i) Mx' > T'p-. The inverse annihilation contribution peaks at T{ ~ 0.28 Mx' and: 
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where y is defined below (51). To a good approximation, y « 292 if > T{ > T'p. On the other 
hand, if {Fix' is very large), y will become suppressed by a factor of exp(—2Mx'/i^max)- 


ii) Mx' < Tp''- The inverse annihilation contribution peaks at Tl ~ T^/1.75, and: 
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Reducing to a Single Sector 


Though the results derived in this Section assume a two-sector cosmology as described in Section it is 
straightforward to reduce these expressions to the single sector case. To see this, we dehne a temperature 
TO = (^-3/2 - yl-4)"^/^ where: 


^0 

max 


( 1 ^/^ / ff.(rRg)V2 

VSy ^ ff*(Tmax) 


1/4 

(MpiR,r2^)V4. 


(59) 


j'O corresponds to the temperature for a given value of A in single sector cosmologies (see eq. (15) in (l^). 
Reducing our expressions to the single sector case amounts to replacing both T and T' with in the above 
expressions for Comparing (59) to (0 -(|20[), this amounts to making the replacements (1 — ?/) —)• 1 

and ri/g'^{T') —)• l/g^{T) in the above expressions. 


4 Implications for UV-motivated Supersymmetric Theories 

In this section, we examine the implications of the results obtained in Section for UV-motivated super- 
symmetric theories that contain moduli fields, and identify regions of parameter space which yield suitable 
DM candidates. As discussed in Section the DM relic abundance in these models is fixed by the following 
parameters: 

Trh, m^,Btot, V, g*{T), g'(T'),Mx', (cru)' (60) 

To simplify our analysis, we will henceforth assume that g*{T) and g'^,{T') are constant, and take g*{T) = 
g'^{T') = 10.75. We also hx g = 0.1, which is a reasonable value assuming the modulus couplings are not 
sequestered from the dark sectoi]^ Relaxing these assumptions will change the computed relic abundance as 
per the formulae in Section |3.2[ but will not qualitatively effect the results presented here. 

As mentioned in Section]^ the parameters Trh, and Btot can be viewed as inputs from the UV theory, 
and are fixed by the couplings and masses of the moduli fields. For a particular UV framework, these quantities 
are constrained to lie within a particular range of values. We will focus here on UV completions which contain 
gravitationally coupled moduli fields while also yielding TeV scale supersymmetry. If the modulus interacts 
gravitationally, dimensional analysis suggests that T^ = and Trh as defined in (j^ is related to 

the modulus mass via: 

Tr„ ^ 14MeV X (61) 

Thus the BBN bound Trh ^ MeV places a lower bound on in the tens of TeV range. 

The range of values for is further restricted by imposing the requirement of TeV scale supersymmetry. 
We focus here on models in which SUSY breaking is mediated to the visible sector via gravitational interactions; 
this arises naturally in theories containing moduli. In the minimal case (i.e. no sequestering or large volume 
suppression of SUSY breaking), the lightest modulus mass is order the gravitino mass m^j 2 , which sets the 
scale of the SUSY breaking parameters For many such models, the scalar superpartner masses will 

be comparable to m^i 2 , while the gauginos may be parametrically lighter by roughly a loop factor. The 
lightest superpartners in the visible sector will then be gauginos whose masses are suppressed with respect 

®For this value of rj, the latest CMB bound on Nf,ff requires that all dark radiation particles have masses greater than ~ 1 eV. 
Otherwise, rj must be smaller. The qualitative features of our results will be the same for smaller g as well. 
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to m^j 2 - This is true for Type II and heterotic models with KKLT-type moduli stabilization, M-theory 
compactifications with stabilized moduli, and also for spectra with pure anomaly mediation. Thus for these 
SUSY models, the requirement of TeV scale supersymmetry along with constraints from BBN imply: 


30TeV<m<^< 0(100) TeV, 


5MeY<TRH<O{l00) MeV, 


(62) 


assuming ci ~ 0(1). This jnstifies our choice of benchmark parameters in ([^. If the DM is an MSSM particle 


there is a tension between (62) and indirect detection constraints, which require Trh > 1 GeV 


The quantity Btot is more difficult to constrain from a theoretical point of view, as it depends on the precise 
interactions between the modnlus and visible/dark sector particles. Nevertheless, if the canonically normalized 
lightest modnlns contains a non-trivial fraction of the modulus that determines the gauge coupling of the visible 
and/or dark sector, then one expects a contribution to Btot by operators of the form f d?6 Wa IU“ where Wa 
is the chiral gauge superheld of either the visible or dark sectoirl Therefore, in M-theory compactifications 11 


and also roughly isotropic Type II compactifications, Btot is expected to be 0(0.1). However, in anisotropic 
compactifications in which the visible and dark sectors are localized at different regions of the internal manifold, 


it is possible that Btot is suppressed, see 17,18 for example. We will consider below a wide range of values 


for Btot to perform as general an analysis as possible. 

In the following, we fix Trh, rrifj, and Btot to particular values, and scan over {av)' and Mx' to give a fairly 
model-independent characterization of the viable regions of DM parameter space. All other parameters are 
taken to their benchmark values ([^. In Figure]^ we have scanned over the {crv)', Mx' parameter space for 
various values of Btot, with Trh = 10 MeV, = 50 TeV for the left column and Trh = 100 MeV, = 150 
TeV for the right column (consistent with the Trh oh scaling manifest in (61)). 

Btot determines both the cross section reqnired for QSEnr (see (34)), and the size of the modulns decay 


contribution (42) in the inefficient annihilation region. Thus the available parameter regions are quite sensitive 


to orders of magnitude changes in Btot- For Btot = 0.1, the viable parameter space effectively splits into two 
regions. In the upper region (av)' > 10“®GeV“^, the relic DM abundance is produced via either QSE^r 
or while in the lower region {(tv)' < 10“^^ and the relic DM abundance is populated via inverse 

annihilations and/or modulus decay. In the inefficient annihilation regime, most of the parameter space with 
Mx' > Trh results in an overabundance of DM due to the modulus decay contribution for Btot = 0.1 (see 


(42)). The value of Mx' where the modulus decay contribution (42) saturates = 0.12 scales like 


Btot~^', thus for smaller values of Btot, much more of the Mx' > Trh parameter space becomes available. 
Particularly, for Btot ^ 10“^ both the FO^°'^ and lAnr mechanisms can give the correct relic abundance for 
a significant portion of the parameter space. These mechanisms are absent for Btot = 0.1, as the DM masses 
required would result in too large a contribution from modulus decay. 


5 Experimental/Observational Consequences 

In this section, we discuss potential experimental probes of the framework analyzed above. As one can imagine, 
since our analysis covers a large range of values for the mass and couplings of DM in Mx' and {(tv)\ there are 
variety of interesting possibilities for observations. A detailed analysis of the various experimental signatures 
which can arise in this framework is beyond the scope of this paper. Instead, we will limit ourselves here to 
making some general and preliminary remarks which will be relevant for futnre studies. 

®This allows the lightest modulus to decay to visible or dark sector gauginos, which would then cascade decay to the DM X'. 
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logic Mx'(GeV) logic Mx'(GeV) 


Figure 6: Left column: scan of the {(tv)', Mx' parameter space with Tb.h = 10 MeV, = 50 TeV, and various values of 
Btot- Right column: similar plots with Trh = 100 MeV, = 150 TeV. All other parameters are fixed to the benchmark values 
1^. Solid (dashed) contours correspond to = 0.12(0.012). Green, blue and gray regions represent < 0.012, 

0.012 < QoAih^ < 0.12 and 0.12 < Qouh^- For these plots we have taken Hi = lO^^F,^, corresponding to ~ 3 TeV (18l. 
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A nice schematic illustration of the framework studied here is provided in Figure From there we see that 
there are three different kinds of couplings, denoted as: {Ay_y, Xy-D and Xd-d}- Now, the very assumption 
that the visible and dark sectors are ‘separate’ sectors implies that the “portal” couplings of type Xy-D are 
parametrically smaller than the {Xy^y, Xd^d} couplings. When this is true, (av) dominantly depends on 
Xy-y, while (av)' depends mostly on Xd-d- However, within this framework, it is the portal couplings of 
type Xy-D that determine the signals for all “standard” searches for dark matter, such as direct-detection, 
indirect-detection, and collider experiments. The portal couplings Xy-D can cover a huge range. At one 
extreme, it is possible to have Xy-D — Xgrav, the latter corresponding to gravitational strength couplings 
suppressed by the Planck scale. In this case, the decay width of the LOSP X, Fx, is comparable to that 
of the modulus P,^ In our work, we have not focused on this case for both theoretical and experimental 
reasons, see Appendix [B| At the other extreme, it is possible that Xy-D is large enough so that the two sectors 
are in thermal equilibrium with each other and thus combine to form one sector. As mentioned above, we 
have also not focused on such a regime. 


Nevertheless, the models considered here can still accommodate a huge range of values 1 Xy-D > Xgrav 
for the portal coupling Xy-D, which in turn allows for a wide variety of DM signals (or lack thereof) in direct 
detection, indirect detection and collider searches. For this range of portal couplings, our results from Section 


3.2 show that the relic abundance for the dark matter X' does not depend on the properties of the LOSP X 


- {Mx,^x, (cu)}, or equivalently the portal couplings Xy-D- Thus, in order to characterize the “standard” 
DM signals which arise in this framework, one must consider explicit dark sector models in which the size 
of the portal couplings Xy-D are calculable. We save this exercise for future work, except for making some 
comments about the consequences of a decaying LOSP. 


The LOSP X, being a visible sector particle, can be produced at colliders. Since it is unstable, it is possible 
that the LOSP is charged and/or colored. Prospects for detecting a charged/colored LOSP at the LHC are 
much better than that for a neutral LOSP, as a charged/colored LOSP will interact with detector materials 
and slow down considerably relative to a neutral LOSP. Charged/colored LOSP decay widths in the range: 
10“^^ GeV > Fx > 10“^^ GeV can be measured in principle. However, subject to model-dependent details, 
large windows in the above range are now disfavored 19 . On the other hand, only decay widths larger than 
around 10“^^ GeV (rx ^ 10“® s) can be measured for a neutral LOSP because then a sizable fraction of 
LOSP particles decay inside the detector. The LOSP decay width Px can be parameterized as: 


rx~^^Mx. (63) 

167r ^ ’ 

Thus, one requires Xy-D ^ 10“® in order for a weak scale neutral LOSP to significantly decay inside the 
detector so that its decay products could be measured in principle. Otherwise the neutral LOSP is stable for 
collider purposes, and manifests itself as missing energy. 

We now describe some possible signatures of the framework that do not depend on portal couplings between 
the visible and dark sectors. 


5.1 Cosmological/Astrophysical Effects 

Here, we comment on astrophysical and cosmological effects arising from two sources - i) that from DM 
couplings of type Xd-d, i-C- from interactions within the dark sector, and ii) from the presence of a modulus- 

^'^This is because the modulus also couples with gravitational strength to both the visible and dark sectors 
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dominated phase in the early Universe. Since these effects are independent of Xv-D couplings, the observables 
which arise are independent of the pattern of “standard” signals for DM. As such, they provide additional 
observables to probe DM and its properties. Some interesting examples of such effects include: 

• Observables sensitive to power spectrum of density fluctuations of dark matter. 

• Observables sensitive to the morphology of galactic DM halos. 

Understanding these and other observables is becoming increasingly important, both because of the realization 
that interactions in dark sector can affect these observables, as well as from the fact that the quantity and 
quality of cosmological and astrophysical data has been getting steadily better. Here, we brieffy discuss the 
following issues: 


Sensitivity to Modulus-Dominated Era: 

The presence of a modulus-dominated era in the early Universe can have important implications. As 
pointed out in 10,20 , this can lead to substantial linear growth of sub-horizon DM perturbations 


during the modulus-dominated era. More precisely, the presence of a (low) reheat temperature sets a 
new cosmological length scale, Lrh = io,RH Hrh)~^, the comoving horizon at the time of reheating. 
Therefore, in the absence of other effects, DM perturbations on length scales I < Lrh grow linearly 
during the modulus dominated phase, and could have interesting observable effects. However, presence 
of other relevant scales can affect whether such sub-horizon growth of DM perturbations are observable 
or not. These scales are described below. 


• Damping of DM Perturbations due to Acoustic Oscillations &: Free-Streaming: 

It is well known that chemical equilibrium is in general different from kinetic equilibrium. In the context 
of DM interactions, the former is set by number-changing interactions in which DM number is not 
preserved, while the latter is set by number-preserving interactions in which DM number is conserved. 
For example, within the standard WIMP paradigm, chemical decoupling leaving to thermal freezeout 
happens much earlier than kinetic decoupling since the interaction rate for the latter is enhanced by the 
relativistic abundance of light SM species in interactions of the type: DM -\- SM —?■ DM SM. 


There are two important scales related to kinetic decoupling that determine the length scale at which 
DM perturbations get damped or suppressed: 


i) Scale arising due to the coupling of DM to the dark radiation fluid (and also to the visible radiation 
and baryons in general). The effect of coupling of DM to visible baryons and radiation is also present for 
standard WIMPs in general 21 , but qualitatively different effects may arise here due to the presence 

It is expected that the DM-DR interactions will give rise 


of dark radiation (DR) in addition p2, 23 


to damped oscillatory features in the DM power spectrum with a characteristic length scale denoted as 
Ld, given by. 


Ld = 


Vkd 
Xd ’ 


(64) 


where pkd is the conformal time at kinetic decoupling, and Xd is a numerical factor of 0(1) (we take 
Xd ~ 7, see for example). 

ii) Scale arising due to the free-streaming of particles after kinetic decoupling. This length scale is 
defined as Ljg = v/adt, where v is the average DM velocity, a is the scale factor, to is the current 
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age of the Universe, and U is a characteristic time which is different for different mechanisms and will 
be discussed shortly. If the universe is radiation dominated at U, then Lfg is given by (see e.g. pO]): 


L 


rad 

fs 


Ja, 




(65) 


where a^q ~ 2.9 x 10 and Hq k, 1.5 x 10 GeV. If the universe is modulus dominated at G, Lfg is 
instead given by: 
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mod 


1/2 

^RH 


HoV^ 




( Mx'Cl 



da + 


a.RH-1 P* —^ Prh 


( 66 ) 


where we have taken H = Hrh {orh during modulus domination. Here a^H corresponds to the 
scale factor at which H = T^, normalized such that a = 1 today. 

Both scales above are present in general, and the damping scale is determined by Lent = max(Lrf, Lfg). The 
scale Lent is relevant in determining the mass of the smallest DM proto-halos: Mproto oc 

As discussed above, DM perturbations on length scales I such that Lent < I < Lrh grow linearly during 
modulus domination and the growth during this era is not washed out by free-streaming and/or acoustic 
damping effects. Thus, these perturbations could have interesting and novel effects. For example, as pointed 


out in 10 , a low reheat temperature of order 10 MeV or so can give rise to an abundance of earth-mass dark 


matter microhalos in the early Universe containing a significant fraction of dark matter. A possible way to 
observe these microhalos is via their strong gravitational lensing effects on quasars |24[|25| , or via their impact 
on pulse arrival times from millisecond pulsars [^. Furthermore, if the portal couplings Xy-D are large 
enough, these DM microhalos can annihilate to 7 -rays, thereby acting as 7 -ray point sources and contributing 
to the 7 -ray background (27-29 . It is worthwhile to explore these possibilities in more detail. 


In the case where Lent > Lrh such that the growth of DM perturbations during modulus domination is 
washed out, the damping of DM perturbations below the scale Lent can still give rise to observable effects. 
Notably, “warm” dark matter with Lent = ~ 1 — 100 Kpe can reconcile many of the discrepancies between 


ACDM cosmology and observations on galactic/sub-galactic scales 30 -32 . If the damping scale becomes too 


large i.e. Lent ^ 1 Mpe, bounds from Lyman-a will start to apply 0 


5.2 Prospects for the Framework 


What can be said about the effects mentioned above vis-a-vis the framework considered? Qualitatively, there 
are two different scenarios which are determined by whether or not is larger than T/,. If < T/,, then 
X' kinetically decouples during radiation domination after the modulus has decayed. Depending on the mass 
and kinetic decoupling temperature, either Ljg or will determine the damping scale Lent- Alternatively, if 
T^d > '^D then the DM kinetically decouples during the modulus dominated phase; Ljg will then determine 
Lent for most of the relevant parameter space. 


In order to discuss observational signatures for the framework considered here, it is pertinent to consider 


what range of values for is expected, given the DM production mechanisms discussed in Section 3.2 
Generically, one expects that crossing symmetry relates the X'X' —)• R'B! annihilation cross section {{av)') 
to the [X'R' — 7 - X'R') elastic scattering cross section (u/J. In the case of fermionic DM annihilating into 
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fermionic B! through a massive bosonic mediator, {<Jv)' (Mx'^ + r'^)/A^ and ~ T'^/A^ where A is the 

mediator mass scale (see e.g. [^[^). If this is the only X' — R' scattering process, X' kinetically decouples 
when the scattering rate Fez ~ drops below the Hubble rate, where u = 1 {T'/Mx') for relativistic 

(non-relativistic) X'. Taking the benchmark parameters ([^ for this example, implies A < 800 GeV 

(250Mx'~^'^^ GeV) if X' kinetically decouples while relativistic (non-relativistic). 

However, in a more realistic model there may be other (e.g. inelastic) processes which also keep X’ in 
kinetic equilibrinm; thus the precise relationship between {(tv)' and is fairly model-dependent. In the 
following, we will treat as a free parameter, though it will be useful to keep the above toy example in 
mind as a benchmark scenario. 


Kinetic Decoupling During Radiation Domination: In this case, the DM particle X' is in kinetic eqni- 
librinm until T' < T'j^. The kinetic decoupling temperature will then determine the length scales in 
(64) and Lfg in (65). Specifically, Lfg is computed using (65) with T'{t^) = and p* = 
where cj = l{Mx'/T^^) for Mx' < T'f^^ {Mx' > F^^om Figure]^ we see that for larger DM masses 

10“^ GeV < Mx' < 10^ GeV and smaller kinetic decoupling temperatures < 0.1 T^, Ljj is larger than 
Lfg and determines Lent and the mass of the smallest proto-halos. All of this parameter space is consistent 
with the upper bounds arising from the observables studied in [^. In the complementary parameter space. 
Lent is determined hy Lfg. A large region of this parameter space is consistent with the Lyman-a forest upper 
bound on Lfg of about 1 Mpc. Finally, for most of the parameter space Lent = iaiax(L^, Lfg) is greater than 
Lrh, implying that growth of DM perturbations in the modulus-dominated era is washed out. Only in a very 
small region of parameter space with 1 < Mx' ^100 GeV and 0.1 < T'j^, one has Lent < Lrh so 

that the memory of growth of DM perturbations on length scales I with Lent < I < Lrh, is retained. This 
can have interesting implications as mentioned previously. 

Kinetic Decoupling During Modulus Domination: In this case, X' kinetically decouples before the begin¬ 
ning of radiation domination such that If kinetic decoupling occurs after X' production, Lfg is 

given by ( 66 ) with T'{t^) = and p* = Note that this scenario requires Mx' > as for 

Mx' < Tf) DM production occurs predominantly when T' < (see Section 3.5). For the allowed parameter 
regions depicted in Fignrej^ one finds in this case that L^ < Lrh < Lfg ^ 1 Mpc, assuming a single DM 
particle accounts for all of the dark matter. 

If kinetic decoupling ocenrs before X' production, T'{a^) is the characteristic temperature at which X' 
production occurs, and p* depends on the mechanism for X' prodnetion. The production mechanisms which 
allow for X' to be produced out of kinetic equilibrium are (see Section 3.2): 


Inverse annihilation: As discussed in Section 3.4.2, DM production from inverse annihilations peaks at 
Ti ~ T^/1.75 for lAr and Tl rs 0.28Mx' for lAnr- We then take p* « \/3T)^ in computing Lfg. 

Production from Modulus Decay: If the DM abnndance comes predominantly from modnlus decay (i.e. 
~ Lldecayh"^ (42)), then and p* « m0/2 assuming 2-body modulus decays. 


To be precise, if X' kinetically decouples before X' production occurs, one mnst replace X' in the Boltzmann 
equations with an integral over the X' phase space distribution function. However for the inverse annihilation 
and modulus decay production mechanisms, terms involving X' can be neglected in the X' Boltzmann equa¬ 
tion; thus the results in Section 3.4 are still valid despite this departure from kinetic eqnilibrinrrp] Nonetheless, 

^^One additional subtlety is that if X' is out of kinetic equilibrium, we are no longer justified in assuming Ex' ~ %/ 3 p '2 _|_ yfx'^ 
in |7|. However if Bx' < 0.1, this subtlety will not significantly effect our results. 
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Figure 7 : Hierarchies among the cosmological length scales Lf^, Ld, Lrh shown in the Mx' — T^d plane, assuming Tl.d < T'd- 
The pink region corresponds to Ld > Lfs, the green region corresponds to Lfs > Ld, the blue region corresponds to Lhh > Lent, 
and the brown region corresponds to Lent > 1 Mpc. The other relevant parameters are set to their benchmark values, see 


in order to properly compute p* and L/s, a precise knowledge of the DM phase space distribution function at 
is required. In lieu of a more precise computation we will use the approximate values for p* quoted above, 
with the understanding that our results for Lfg are meant to be qualitative. 

Figure [^summarizes the cosmological length scales which can arise in the case where X' is produced out 
of kinetic equilibrium. Because X' is not coupled to the dark radiation bath when produced, there is no 
acoustic damping effect to consider; thus Lent = We see from Figure [^ that for the I Ay case, most of the 
parameter space easily avoids Lyman-a constraints. The lA^ scenario can also naturally accomodate warm 
DM candidates, with Ljs ~ 1 — 100 kpc. Perhaps more interestingly, we see that for a majority of the lAnr 
parameter space, Lfs < Lrh- Thus, the linear growth of DM perturbations during modulus domination is 
not washed out for a large portion of the lAnr parameter space, leading to potentially interesting effects as 
discussed above. 


Finally, let us comment on the case where relic DM is produced from modulus decay. If DM particles in 
this scenario are kinetically decoupled at T'jj , they will be highly boosted when produced from modulus decay. 
If the modulus decay contribution is the dominant contribution to the overall DM abundance, DM masses 
within the range 10“^ GeV < Mx' ^ few GeV are at odds with Lyman-a bounds < 1 Mpc; this is evident 
from Figure Thus if the relic DM is predominantly produced via modulus decays, Lyman-a constraints 
require Mx' > 0(1) GeV; this in turn implies Btot ^ 10“^ as can be seen from (42). 


To summarize, we find that there are various interesting possibilities for cosmological/astrophysical ob¬ 
servables which can probe the framework considered, both in terms of providing constraints on the parameter 
space as well as by providing insights for potential signals. In particular, we find that there is sensitivity to 
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Figure 8: Plot of the free-streaming length L/s in Mpc. The blue lines correspond to DM produced predominantly via inverse 
annihilations, while the black lines correspond to DM produced predominantly via modulus decay. The solid lines were obtained 
for Trh = 10 MeV, = 50 TeV while the dashed lines were obtained for Trh = 100 MeV, = 150 TeV. The other relevant 
parameters are chosen as in 


the modulus domination era for a large portion of the lA^r parameter space. This is in contrast to the result 
obtained in 20 , primarily because the framework considered here encompasses a wider variety of DM masses 
and couplings compared to the analysis in pO 36 . The results obtained in this section are largely qualitative. 


It would therefore be interesting to carry out a more detailed and comprehensive analysis of the constraints 
and potential observations which have been suggested in this section. 


6 Hidden/Extra Sectors in Explicit String Constructions 

The system of Boltzmann equations studied in this paper are relevant for cases in which the the dark matter is 
located in a hidden sector, weakly coupled to the observable sector, and the universe undergoes a substantial 
period in which the total energy density is controlled by a single modulus field. In addition, we have chosen 
to study a benchmark case in which the number of relativistic degrees of freedom for the dark sector is similar 
to that of the Standard Model. How well motivated is this framework? 

Certainly, the presence of hidden sectors is generic in string theory, as observed in the earliest days of 
model-building based on the Eg x heterotic string. Furthermore, every four-dimensional effective supergrav¬ 
ity theory representing a string compactification has moduli. Barring some remarkable feat of engineering, 
therefore, a long period of modulus domination in the early universe is essentially guaranteed. These two 
components of our framework are therefore exceedingly well motivated from the point of view of string theory. 

But this paper adds a third component: a portal that connects the hidden sector and observable sector. 
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Our analysis generally assumes that the strength of this presumed coupling is greater than that of gravity- 
mediated operators suppressed by the Planck scale. Interactions between the dark sector and the observable 
sector complicate the system of Boltzmann equations, as we have explained at length. Measurements neces¬ 
sarily constrain the observable sector, but the presence of a portal of appreciable strength mean that these 
measurements also constrain the nature of the dark, hidden sector in a manner that would not exist if the 
coupling between the sectors was utterly negligible. This constraint is summarized in the very first equation 
presented in the paper, and it involves the number of relativistic degrees of freedom in the hidden sector. It 
is thus important to ask what, if any, statements can be made about the nature of hidden sectors in actual 
string constructions, and what sorts of interactions are observed to exist between these hidden sectors and 
the Standard Model. 


6.1 Heterotic Orbifolds 


Calculating the massless spectrum in a string compactihcation is easiest to perform in cases where conformal 
field theory tools are available. This tends to restrict explicit calculations to orbifolds and their orientifold 
analogues. These techniques have been used extensively in weakly-coupled heterotic string theory, but also in 
Type II string theory. The latter have generally been conducted in the context of Type IIA theory compactified 
on orientifolds with intersecting Dg-branes. 


String phenomenologists tend to be concerned primarily with the observable sector, and ensuring that 
three generations of fields charged under the Standard Model gauge group emerge from the compactihcation. 
The hidden sector is often left undetermined, or only computed years later when they become necessary (for 
example) in guaranteeing global consistency conditions. Therefore, meaningful examples in the literature 
are relatively sparse. An important early computation involved Zq asymmetric orbifolds of heterotic string 
theory, in which quasi-realistic GUT models were constructed at higher Kac-Moody level 37,38 . A search for 


constructions which yielded an Eq, 50(10), SU{6) or SU{5) GUT model was conducted. Satisfactory cases 
were required to have three (net) families of GUT representations capable of realizing the SM, and an adjoint 
Higgs representation for breaking the group to the Standard Model. The three families in this case arise from 
demanding a outer-automorphism. 


In addition to the GUT gauge group, the hidden sector groups were identified, and the massless matter 
content for all sectors was computed. For Eq and 50(10) GUTs, the hidden sector consisted of at least one, 
and sometimes two, SU{2) factors. In one class of constructions bifundamental representations between the 
SM and a hidden 517(2) were identified, suggesting the possibility of a Higgs-like portal between the sectors. 
In addition, there were several (non-anomalous) 17(1) factors and the SM states typically carried charges 
under these ‘hidden’ 17(1)’s. For the 51/(5) and 51/(6) GUTs, the hidden sector gauge groups can again be 
51/(2), but occasionally 51/(3) and even 51/(4) were observed. Again, some twisted sectors tend to contain 
states that are fundamentals of the SU{N) GUT group, but also a doublet under a hidden 51/(2) factor. We 
note that in most of the cases studied, the 51/(2) factor remains weakly coupled to very low energies. This 
suggests that a relatively sparse dark sector that contains a WIMP, interacting with the observable sector via 
a Higgs portal or 1/(1) portal, would not be unusual in this particular class of theories. 


An example in which the hidden sector analysis followed that of the observable sector by nearly a decade, 
is the venerable case of heterotic Eg, x Eg string theory compactified on a Zg orbifold. The original work of 
classifying all possible twist embeddings for the Zg orbifold, with two Wilson lines, that yield the Standard 
Model gauge group was performed in the late 1980’s 39 . But this analysis did not fully specify the embedding 
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for the hidden sector factor. This completion was performed in 2000 40 . The classification required three 


generations of Standard Model matter. In practice it is the demand that three quark doublets be present that 
puts the most restrictions on the allowed Wilson lines. This, in turn, restricts the allowed hidden sector gauge 
group to the relatively small list of 50(10), 51/(5), 51/(4), 5t/(3) and SU{2), plus additional U{1) factors to 
fill the rank-eight product group (one of which will be anomalous). In a follow-up study [^, the field content 
charged under the full rank-sixteen gauge group was computed, and it was common to find states in twisted 
sectors which were bifundamental between the SM SU{2)l and a hidden 51/(2). In addition. Standard Model 
fields and hidden sector fields were generally charged under any number of common 1/(1) factors. 


Years later, a much more exhaustive search was performed, this time in the context of the Zq-u orbifold 
of heterotic string theory. This so-called ‘mini-landscape’ study |42| required a gauge embedding of the 
orbifold action such that an intermediary 50(10) or Eq GUT structure emerges. The authors then scan over 
all possible completions of the embedding with up to two Wilson lines such that intermediate GUT gauge 
group ultimately breaks to the SM gauge group. Further requirements included the demand of three (net) 
generations of SM fields and a hypercharge candidate which is non-anomalous. Unfortunately, only a single 
explicit example of a hidden sector was given, that of an 50(8) x 51/(2) hidden sector. Again, the massless 
spectrum contained states which were bifundamental under the SU{2)l of the Standard Model and a hidden 
SU (2), though these states were vector-like with respect to the overall gauge group and might therefore receive 
large masses if an appropriate set of singlet vacuum expectation values were to arise. A follow-up study relaxed 
the restriction to intermediate 50(10) and Eq structures, and allowed up to three Wilson lines [^. Once 
again, however, a single example of a hidden sector was illustrated, containing both an 51/(3) and an 5t/(5) 
factor. Interestingly, this example had states which were fundamentals under the hidden SU{3) and yet were 
charged under various parts of the Standard Model gauge group. It is hard to assess just how generic such 
portals are in this promising class of constructions, as the raw data was not presented in the papers. However, 
a later paper by Goodsell et al. investigated this same data set, analyzing the prospects for kinetic mixing 
between U{1)y and hidden sector t/(l) factors. The authors found that over 95% of the models allow for 
such mixing, and some explicit dark sectors and dark forces were constructed. Their conclusion was that such 
sectors and portals were indeed ‘generic’ in this class. 


It would be interesting to know if such properties were also common in smooth Calabi-Yau compacti- 
fications of the heterotic string, away from the orbifold point in moduli space. A systematic investigation 
of hidden sectors in this context has yet to be performed. An initial foray into the subject was presented 
in [45] , in which a search was conducted for consistent vector bundle configurations of the Eg x Eg hidden 
sector, given a holomorphic observable sector bundle with structure group 5t/(4), which was shown to allow 
for the three-generation Standard Model held content [^. From this, two examples were presented, with 
gauge groups 50(12) and Ej, and neither case seemed to contain a portal between the two sectors. But we 
note that this paper was meant as a proof-of-concept, not an exhaustive survey. 


6.2 Local Models in a Global Embedding 

Recent years has seen an explosion in model building in the context of Type II string theory compactihed on 
orientifolds. For the most part, this model building has occurred in the form of ‘local’ models: the study of 
ZZ-branes at singularities in the Galabi-Yau manifold, in which other effects (including the presence of possible 
hidden sectors) can be safely neglected. Models are constructed using representative quiver gauge theories, 
or the related techniques of dimer diagrams/brane tilings. 
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Unfortunately, these studies tend to focus exclusively on the Standard Model field content. Hidden sec¬ 
tors only emerge when an effort is made to embed these local constructions in a global Calabi-Yau context. 
To a hrst approximation, such global embeddings amount to the imposition of certain consistency condi¬ 
tions, including the requirement of iV = 1 supersymmetry, Ramond-Ramond tadpole cancellation and various 
anomaly constraints. These additional requirements generally necessitate sectors beyond the local Standard 
Model quiver (i.e. hidden sectors), as the Standard Model theory generally does not satisfy them on its own. 
In fact, in the context of quiver gauge theories, these consistency conditions will generally require that matter 
charged under hidden sector ‘nodes’ are binfundamental with the nodes of the SM gauge group. 

An early example involved Type IIA orientifolds on T®/Z 2 x Z 2 with intersectiong Hg-branes 47 . These 
models achieve the three-generation Standard Model via the Pati-Salam gauge group S'f7(4) x SU{2) x SU{2). 
A scan was performed over all possible brane configurations and wrapping numbers consistent with the 
Standard Model field content (via the Pati-Salam symmetry) and global consistency conditions. This yielded 
explicit hidden sectors which could then be classified. Typical hidden sectors involved the USp{4) and USp{2) 
symplectic groups. In some cases, the field content charged under these gauge groups allowed for confinement 
of the USp{n) gauge group, and the authors speculate as to the appearance of various ‘mesonic’ and ‘baryonic’ 
bound states. It is noteworthy that such composites would generally carry charges under the various residual 
[/(!) symmetries, including that of the Standard Model. Typically, the number of such objects in the massless 
spectrum was of order ten, consistent with the number of degrees of freedom in the Standard Model below 
the QCD confinement scale. 

A more expansive survey was conducted some time later 


48 . In this case the survey began with the 


original three-node ‘Madrid’ quiver 49 and all its three-node generalizations. These quivers represent gauge 
theories which contain the field content of the Standard Model. From this, additional nodes (i.e. gauge groups) 
were added to the quiver until all global embedding conditions were satisfied. All of the quiver extensions 
considered in this paper had an anomalous 17(1) factor, under which the newly introduced ‘hidden’ states are 
chiral, not vector-like. As a result, there is a mixed anomaly between this 17(1) factor and any other 17(1) 
under which these states are charged. In particular, when the hidden sector states carry hypercharge, a mixed 
anomaly between hypercharge and the anomalous 17(1) provides a portal between the two sectors, with a 


potentially light mediating Z' boson. For the phenomenology, and some toy models, see 50 


In Type IIB string theory, much of the recent work has focused on singularities of toric del Pezzo surfaces 
(dPn surfaces). The advantage here is that toric surfaces afford a certain ‘modularity’ in constructing models. 


in which one can work in a bottom-up approach, beginning with various phenomenological demands 51 
Another benefit is access to the large Kreuzer-Skarke database of reflexive polytopes 52 , which generate 


these toric ambient spaces, and the Calabi-Yau manifold realized as a hypersurface within these ambient 
spaces (^ . 

Local models with promising phenomenological features were constructed in this context in recent years |54l 


55 . While these early efforts concentrated almost exclusively on the observable (Standard Model) sector, some 


attempts at embedding these into a global Calabi-Yau context have been made, by enforcing consistency 


conditions such as the Ramond-Ramond tadpole conditions and vanishing of K-theory torsion charges 56,57 


In these papers, some rudimentary hidden sectors were constructed. The authors chose flux parameters in 
such a way as to avoid chiral matter in the hidden sector, resulting in a pure 50(8) or 5C/(4) gauge theory 
with no portal to the observable sector. However, this was again a proof-of-principle and not an exhaustive 
scan over all possible hidden sector conhgurations. 

To our knowledge, no such survey has been conducted within the Type IIB context in analogy with 
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the above-mentioned work in Type IIA. However, some interesting examples of non-trivial hidden, or ‘dark’ 
sectors, were constructed using the ‘toric Lego’ approach of [51] . One such example involved the construction 
of a “dark sector” which mimics the MSSM (visible) sector. This was a toy model designed to exhibit the 
power of the modular approach. The model was based on two dPo singularities and a dPi singularity - the 
former pair for the visible and dark sectors, the last for the SUSY breaking sector. The dark and visible sectors 
were patterned on the phenomenological model of 58 , in which kinetic mixing between U{1) factors in the 
observable and dark sectors provide the portal. The global embedding was identified some time later 
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by identifying those reflexive polytopes from the Kreuzer-Skarke database with the appropriate singularity 
structure in one of their two-dimensional faces to give rise to this trio of sectors. Remarkably, nearly 300,000 
such polytopes were shown to exist, implying at least as many (and perhaps many more) Calabi-Yau manifolds 
which would generate this model upon compactification. 


Analysis of hidden sectors in bona fide string constructions - at the level needed to describe early universe 
dynamics - is still in its earliest stages, lagging the construction of viable observable sectors in many respects. 
The areas that have been investigated were those that were identified as being phenomenologically interesting 
from the point of view of observable sector physics, and may not be a representative sample of heterotic or 
Type II string theory, let alone the entire string theory landscape. Nevertheless, the basic elements that are 
needed for our cosmological framework are often present. 


7 Summary and Future Directions 

In this work, we have provided a general classification of dark matter models in a Universe which undergoes 
a phase of pressure-less matter (modulus) domination. Such non-thermal cosmological histories are predicted 
in a wide class of UV completions to the Standard Model (e.g. compactified string theories), and are also phe¬ 
nomenologically viable provided that the matter dominated phase ends before BBN. Our analysis generalizes 
previous works by going far beyond the standard WIMP paradigm. In particular: 

• We consider DM masses and annihilation cross sections which span several orders of magnitude above 
and below the electroweak scale. 

• We allow the possibility that DM in thermal equilibrium with a ‘dark sector’, whose temperature need 
not be the same as that of the visible sector. 


Upon analyzing the relevant Boltzmann equations, we classify the mechanisms by which relic DM can be 
produced. We find four distinct mechanisms (QS'E'nr,/A|j.^nr} and each of which have 

different parametrics for The first three mechanisms are different from standard thermal freeze-out. 

We derive semi-analytic approximations for these various production mechanisms, and discuss their regimes 


of validity. For the convenience of the reader, these results are summarized in Section 3.5 


Our results have interesting implications for supersymmetric theories containing moduli fields. As discussed 
in Sections]^ anddoes not depend on the masses or couplings of the (unstable) lightest visible sector 
superpartner (LOSP), provided the LOSP decays before the end of modulus dominatior Once the modulus 
mass and couplings are fixed and the dark relativistic degrees of freedom g'fiT') are specified, depends 

only on Mx' and {crv)'. In Section we fixed the modulus mass and couplings by considering models with 


^The contrary case is briefly considered in Appendix p3^ 
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gravity mediated SUSY breaking in which is of order the gravitino mass. We mapped out the parameter 
space of these models by scanning over M^', {c^v)' for various values of Btot, see Figure]^ Here Btot is 
the branching ratio of the modulus decay into DM, including contributions from intermediate states. For 
Btot ~ 0(0.1), the viable DM parameter space splits into two seperate regions: large annihilation cross 
section {av)' > 10 “® GeV“^, or small annihilation cross section {(tv)' < 10 GeV Intermediate values 
of {cFv)' result in DM overproduction. Moreover in the {(Jv)' < 10“^^ GeV“^ region, the DM mass must be 
< 100 MeV to avoid being overproduced by moduli decay. If however the modulus branching ratio to DM is 
suppressed i.e. Btot ^ 1, much more of the DM parameter space becomes available. These features can easily 
be inferred from Figure]^ 

We have also briefly discussed potential experimental signatures for the theoretical framework considered 
here. Since is insensitive to the portal couplings between the visible and dark sectors for the models 

considered, the “standard” DM signals in direct detection, indirect detection and collider experiments, which 
crucially depend on portal couplings between the visible and dark sectors, can cover a wide range of possibilities 
are rather model-dependent. On the other hand, observables which involve couplings within the dark sector 
yield more robust predictions, as these couplings are correlated with the DM relic abundance. One such set 
of observables involves the power spectrum of DM density perturbations. If the DM kinetically decouples 
during the radiation dominated era after BBN, the sensitivity of DM density perturbations to the modulus 
dominated phase is maintained only for a very small region of parameter space, as shown in Figure On the 
other hand, when DM kinetically decouples during modulus domination, the power spectrum of DM density 
perturbations depends on the mechanism by which relic DM is produced: 


If DM is produced by annihilation of thermal bath particles while the DM is non-relativistic (we call 


this case see Section 3.5), the free-streaming length is smaller than the comoving horizon at Tjuj. 


The linear growth of DM density perturbations during modulus domination is not washed out, leading 


to potentially interesting astrophysical signatures as discussed in 10 


If DM is produced by annihilation of thermal bath particles while the DM is relativistic (we call this case 


lAj-, see Section 3.5), the free-streaming length is larger than the comoving horizon at Tjih- Even though 


the growth of DM perturbations during modulus domination is erased, a large region of parameter space 
yields Lfs ~ 1 — 100 Kpc which leads to signatures similar to warm DM. 

If DM is dominantly produced by modulus decay, then the DM has large free streaming lengths Lfs > 1 
Mpc, which is in tension with constraints on warm dark matter from Lyman-a measurements. 


There are many opportunities for future research. From the point of view of the Boltzmann equations, 
including n —>■ 2 annihilation processes where n > 3 would be worth understanding in this framework (this 
would be the non-thermal analog of [60]). From the point of view of model-building, it would be worthwhile 
to study explicit models of DM candidates and portal interactions within the general framework so that 
detailed predictions for “standard” DM signals (e.g. direct and indirect detection) could be made. From 
a string theory perspective, although there is some existing work on dark sectors and portal interactions 
as described in Section clearly much more needs to be done. Finally, our discussion in Section of the 
astrophysical/cosmological effects of DM interactions within its own sector has been largely qualitative. A 
more precise analysis would involve solving for the DM phase space distribution at kinetic decoupling in order 
to determine the appropriate transfer function relevant for the power spectrum of DM density fluctuations. 
We hope that future studies in these directions will help shed important light on the nature of dark matter. 
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A Justifying Approximations for R and R' 


In Section 3.1, analytic approximations for R and R' were obtained assuming that all other terms aside from 
the modulus decay term can be neglected in dR'/dA and dR/dA if Mx',Mx 'C m^. In this appendix, we will 
justify this approximation. Note from Q that the modulus decay terms in dR'/dA and dR/dA grow like 
during the modulus domination phase, and peak when T ^ T^. Thus in determining whether or not certain 
terms in dR'/dA and dR/dA are negligible compared to the modulus decay term, it is sufficient to focus on 
the Boltzmann equations at temperatures near Tq. 

First, consider the X —?• A' + ... decay term in dR'/dA and dR/dA. At T > T^, X has already reached 
QSE, assuming (Tx) > T^. Taking X = Aqse with 6 « 1, the A' —)• A + ... decay terms are given by: 


;dR 

dA 


= BxBx^x'rcV^ 


) 


77 -^ = BxBx^x'R' 


Ex — Ex' 


A3/24>+... (67) 


Here Bx^x'r and Bx^x'R' are the branching fractions of A into X'R and X'R'. Thus we see that the 
A — > A' + ... decay terms are suppressed with respect to the modulus decay term by a factor of {Ex—Ex’)/rn^; 
a similar conclusion holds if A does not decay to X'. Next, consider the annihilation terms. For the {av) term 
in dR/dA, Aeq ~ 0 and A « Aqse for T ~ Tp. Thus for temperature-independent {av), the annihilation 
term in dR/dA falls like for T >Td, and will be numerically insignificant at To due to suppression by 

negative powers of the scale factor. 

The argument for the {av)' term in dR'/dA is less straightforward. First, consider the case where {av)' > 
{av)'^ such that X' reaches QSE at T' ~ T^. If Aqg^ ~ -^eq) annihilation term vanishes and is trivially 
negligible. If instead A'q is negligible in Aqgg (32) at T' ~ (as is the case for QSEnr), we can take 
X' PS X'qg^ and write the {av)' term as: 


TjdB' _ 1/2 n 

dA~^P 


) 


( 68 ) 


which is suppressed with respect to the modulus decay term by a factor of Ex'/rn^. Now consider the case 
where {av)' < {av)'^ such that X' is not in QSE at T^. As discussed in Section 3.4 we can write X' at 


r'>r^as a' = a;„,, + a'. 


where comes from integrating the modulus decay term: 


X' . = -c 1/2 

^mod "'"O 


TrhB tot 


^3/2^1/2 


(69) 


and is determined by the {av)' term. In the case where Mx' > Tr/, A ^„„ is negligible compared to 

X'^^^ unless Rtot ^ 1 (see Sections 3.4.1 and 3.4.2). Taking X' ps A(^q^, the {av)' term in dR!/dA can be 
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written as: 


(70) 


H 


dR 

~dA 


" = i^A3/2c^ + ... 

^ 9 ^ V rricj) J {av)[. 


Thus in the case where {(tvY < {crv)'^ and Mx' > T'j^, the {avY term in dR'/dA is suppressed by at least a 
factor of Ex'/nitf) with respect to the modulus decay term. 

Finally, consider the case where {avY < {(^vYc corresponding to the lA^ scenario (see 

3.4.2). In this case X' <C and we can write the {crvY annihilation term as: 


Section 


H 


JR 

~dA 


1 1 p‘ ' 57r65^r')2r'2 


+ 


(71) 


where we have used to relate A and T'. Evaluating ( |71[ ) at T' = T^, we obtain: 

dR! /^r, 


H 


dA 


... 1/2 ^ ( d 9 *{ TRH )\ ' / 48\/3c^ 2^(3)2 (an)'\ 3/2 , 

. - c. A4, r„« , (^^ j ( 5,«.;(r') j - 

_ 1/2 ( 0-16cgC(3) \ /3/2 


(72) 

(73) 


Thus the {(JvY term in dR'/dA is suppressed with respect to the modulus decay term by a factor of the order 
of 0.1 {(tvY / {o'vYc- 

To summarize, the above arguments show that the approximations made in solving the equations for R 
and R' in solving Q are justified, as can also be confirmed by the agreement of the approximate and exact 
solutions in Appendix 


B A Very Long-lived X Particle (Fx T^) 

For most of this work, we have assumed Tx > 0(1) T^ such that X decays are efficient before the end of 
modulus domination. This assumption is well-motivated from both theoretical and phenomenological points 
of view. To see this, note that the modulus decays through Planck suppressed operators such that the decay 
width is parametrically given by: T^ ~ rsj 10 GeV for rsj 50 TeV. Thus, as long as the visible and 

dark sectors are coupled by larger than gravitational strength interactions, one expects Tx F^ for a wide 
class of dark sector models. This is also true if the coupling between the two sectors arises by integrating out 
Kaluza-Klein (KK) modes of the extra dimensions or heavy GUT multiplets of some underlying GUT model, 
as even these mediators are lighter than the Planck scale. In addition, from a phenomenological point of view, 
X decays to visible sector particles can spoil the successful predictions of BBN ifFx < H{Tbbn) rsj 
where Tbbn rsj 1 MeV [61] . To avoid these constraints, for Fx ^ F^, Fx should lie in a narrow window: 

c. F</Mpi > FxMpi > T^bbn- (74) 

Despite these considerations, for completeness we briefly discuss in this appendix the case where Fx ^ F^. 
In this case, X is effectively stable during modulus domination (as H > F,/ > Fx)- Thus for 77 > Fx we 
can treat X as a stable relic. If X is a WIMP, its comoving abundance will become fixed at T ~ Tb via 
the QSEnv mechanism, which is the precise generalization of the non-thermal WIMP miracle [^[^. Once the 
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Hubble parameter drops below Tx during radiation domination, the remaining X abundance will decay to 
yield X' particles. The dynamics of such a process was studied in detail in [^. From the results of [^, we 
see that there are three possibilities for the resulting parametrics of VLx'h?'- 

• X' is in equilibrium when H = Tx (which is only possible for and FO\^). X' will continue to 

track its equilibrium abundance until freeze-out. In this case Qdm is completely insensitive to X 
decays. 

• X' is out of equilibrium when H = Tx, and X decays yield an X' abundance which is less than the 
critical abundance required for X' annihilations. This gives rise to the freezeout &: decay {FOSzD) 

In terms of dimensionless comoving variables, Xqse(^c) < ^crit|_H'=r ’ 

is defined 


mechanism described in 


16 


where Xqse(^c) is given oy the QSEnr mechanism as described in Section 


3.3.1 


and X;^it 


in (23). The resulting contribution to the X' comoving abundance is insensitive to Tx, and is given 
simply by AX' « Xqse(^c)- This contribution must be added to the X' abundance which results from 


the production mechanisms described in Section 3.2 


X' is out of equilibrium when H = Tx, and X decays yield an X' abundance which exceeds the 
critical abundance required for X' annihilations. In terms of dimensionless comoving variables this 
occurs if Xqse(Xc) > The X' particles produced from X decays will then annihilate until 

X' Ri X'I. This was referred to as the ^^freezeout & decay and re-annihilation" (FOSzO^) in 
the resulting X' relic abundance scales like Tl om oc ^ 


16 


-. 1/2 


{(Tv)' 


m 


Before concluding this appendix, we remark that the ^freeze-in’’ mechanisms (FI and FF) described 
are not important for the models considered here. Recall that FI is due to X —)• X' + ... decays 


16 
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which occur during the radiation domination era when X is still relativistic and in equilibrium. However, it 
turns out that freeze-in due to X decays is negligible during the modulus dominated era. To see this, consider 
the X decay term in dX'/dA. We saw in section 3.1.3 that X attains QSE at some scale factor A (say Ax) 


before A^ if Tx > 0(1) T^. For A < Ax, X is given by X « Xeq, while for Ax < A < Ao, X is given by 


X ~ Xqse. In the analysis in Section 3.1.3, the effect of X decays when 1 < X < Ax and X ~ X^q, which 


corresponds to freeze-in effects from X decays, was neglected. To see that it is justified to do so, note that 
the integration of the decay term gives (up to overall constants): 


r/iD r2lx r^i-D 

I (iXXX^/2_ / dAXeqA^/^+ / dAXQSEA^/^ 

O dAArv + Iao^i^ (bxBx^Fk) 

1 3 V ^xm^J 






(75) 


where Xx corresponds to the scale factor at which either X becomes non-relativistic or X enters QSE 


(whichever occurs first). Comparing the first and second terms in (75), we find: 


dXXeqXV^ ^ fToV ( Txm^ \ ^ /2|^£xm^\ 


/^^dXXQSE ^^/2 \Tx) \BxT^Tl,qj) Tx^^i \ Tx^Bx ) 


(76) 


where we have used T and (^Tmax)®/^/ ~ Trh^ (see (18)). There are now two possibilities 

for Tx- If X enters QSE before X becomes non-relativistic, then Tx ~ (FxMpiTR//^)^/^ > Mx- If instead 
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X becomes non-relativistic before QSE is reached, then Tx ~ Mx and Fx < 


Mi 


the latter case, the ratio (76) is maximized for Tx ~ Mx, and one gets: 




Since Tx is smaller in 


10 


-13 


dA Xqse ^^/2 Bx Ml Bx \ 10 MeV 


Trh 


( ^<t> \ f 100GeV \ 

UOOTeV; \ Mx ) ' ^ 


Thus the freeze-in production of X' from X decays can be neglected for reasonable choices of parameters, 
provided Bx is not extremely tiny. 


C Accuracy of Approximate Solutions 



Figure 9: Left: QoAih^ as a function of {(tv)' for Mx' ~ 10 GeV and Btot = 0.1. The green curve shows the numerical 
solution, while the black curve shows the approximate QSEm solution (34l. The vertical dashed line represents {(tv)' = as 

defined in ( |34[ ), while the horizontal dashed line shows represents the modulus decay contribution given in \42\ , which is valid for 
{av) < Right: the ratio of the approximate result for QSEnr to the exact result. 


In this section, we compare the semi-analytic approximations obtained above with the full numerical 
solution to 0. The accuracy of these approximations is depicted in Figures and In these plots we 
use the benchmark values of parameters as in ([^; however we take g'^,{T') = 20 7 ^ 9*{T) to ensure that the 
g'^{T') dependence has been properly captured. Figure shows the accuracy of the approximate solutions 
for Mx' = 10 GeV and i?tot = 0.1. In the left plot, the green curve shows the numerical solution, while the 
black curve in the left plot shows the approximate expression for DM production through QSEnr (37). The 
right-hand plot shows the ratio of the approximate QSEnr result to the exact result, which is close to unity 
if (av)' {o'v)l 
Figure 


10 


shows the accuracy of the approximate solutions for Btot = 0; note that in this case {crv)'^ 
becomes effectively infinite for Mx' > T^, see (34). The left plot shows as a function of {av)' for 

Mix' = 10 GeV. For these parameters, DM production occurs either via for {av)' > ((Tu)q or via 

lAnr for {av)' < (ctu)q where ((Tu)q is defined in (45). The green curve shows the numerical solution; the red 
curve shows the approximation for lA^r (50); and the black curve shows the approximation for (44). 

The right plot shows a similar plot with Mx' = 10“® GeV. In this case DM production occurs via lA^ for 
{av)' < {av)'^ and via thermal freeze-out {FOl^"^ and for {av)' > {av)'^. The green curve shows the 
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Figure 10: Left: Qouh^ as a function of {cv)' for Mx> ~ 10 GeV and Btot = 0. Right: similar plot for Mx' = IQ-® GeV 
and Btot = 0. The green curves show the numerical solution, the red curve shows the approximation for lAm (left) and lA,- 
(right), while the black curve shows the app roximation for _F'O(0°‘i (left) and (right). In the left plot the vertical dashed 

line represents ( av)' = (o-v)g, defined in (45 I, while in the right plot the vertical line represents (o-v)' = (crv)'^ in the case where 
Mx' < T'j^ (see (341). The dashed horizontal line in the right plot shows the approximate solution for FOl^'^. 


numerical solution, the red curve shows the approximate expression for lA^ (51), while the black curve shows 
the approximate expression for (39). 

are accurate to within 


Within their respective regimes of validity, (50), (51) and (44) 
~ 15%. 


5%, while (40) and (39) are accurate to within 


D Temperature Dependence of {av 


In Section (semi)-analytic expressions for were obtained assuming that (av)' is temperature 

dependent. In this section, we generalize the results of Section 


3.2 


for temperature dependent {crv}'. For 
scenarios where the contribution to Udm is determined by non-relativistic X' annihilation {QSE'^'^, FO^'^, 
and lAar) we will consider p-wave annihilations where {crv)' = T' / Is?. For scenarios where the 
contribution to is determined by relativistic annihilations {IA? we consider the case where {(Jv)' = 

T^/A^, corresponding to annihilation through a heavy bosonic mediator. Note that for is 


independent of {crv)' so (40) holds regardless of the temperature dependence of (crv)'. 


For QSFar, (36) and (37) are still valid for p-wave annihilation, provided the annihilation cross section 


is parameterized as: 


, r / Ad 

=X3 = V X 


(78) 


where in the second equality we have assumed T' ex A ^ as in radiation domination. In order to match 


the numerical result, we instead use k = 1.8 in (36) and (37). 
For FO^^, the expression for x'p is given by: 






(79) 
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while ^Duh'^ is given by: 
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Here = 2 is a constant which is chosen to match the full numerical result. 


For the expression for x'p is given by: 

( ^ 


x'p = In 

while is given by: 
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Here = 5/4 is a constant which is chosen to match the full numerical result (see also 14 ). 


(81) 


(82) 


For I Aar, it is straightforward to show that for {crv)' = T'/M^, the expression analogous to (50) is given 
by: 


f^ann ^ [.f Anr] 
where Xp is given by: 
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(83) 


(84) 


The integrand peaks at T/ « 0.33Mx'', if 7/ T^ax T/ ;i> T'p, Xp ^ 80. 

For lAr, we are interested in the case where {(tv)' = T'’^ j (see above). The expression analgous to (51) 
is given by: 
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48c52c(3)2 
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(85) 


In the above, T/ ~ T'p/1.35 is chosen to match the numerical result, and is related to the temperature at 
which the integrand of f dAT'^A'^^^n'gq^ff~^ peaks. Note that we can recover (85) from (51) by making 
the replacement: 


(crv)' —>• 0.17 X 


[ A^gUT^y/2 J RH 


( 86 ) 
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